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ABSTRACT 

We investigate properties of iron fluorescent line arising as a result of illumination of a 
black hole accretion disc by an X-ray source located above the disc surface. We study 
in details the light-bending model of variability of the line, extending previous work 
on the subject. 

We indicate bending of photon trajectories to the equatorial plane, which is a 
distinct property of the Kerr metric, as the most feasible effect underlying reduced 
variability of the line observed in several objects. A model involving an X-ray source 
with a varying radial distance, located within a few central gravitational radii around 
a rapidly rotating black hole, close to the disc surface, may explain both the elongated 
red wing of the line profile and the complex variability pattern observed in MCG-6- 
30-15 by XMM-Newton. 

We point out also that illumination by radiation which returns to the disc (fol- 
lowing the previous reflection) contributes significantly to formation of the line profile 
in some cases. As a result of this effect, the line profile always has a pronounced blue 
peak (which is not observed in the deep minimum state in MCG-6-30-15), unless the 
reflecting material is absent within the innermost 2-3 gravitational radii. 
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1 INTRODUCTION 

Broad iron lines observed from many black hole systems 
most likely originate from the innermost regions of an accre- 
tion disc and their profiles are shaped by gravitational red- 
shift and Doppler shifts. Modelling of the lines observed in 
several objects requires strongly enhanced fluorescent emis- 
sion from a few gravitational radii (e.g., Wilms et al. 2001, 
Fabian et al. 2002, Miller et al. 2002, Miller et al. 2004; see 
review in Reynolds & Nowak 2003), which in turn indicates 
that a primary source of hard X-ray emission must also be 
located close to the black hole. Thus, both the primary and 
reflected emission should be subject to strong gravity ef- 
fects. These effects are also tentatively considered as an 
explanation of the complex variability pattern characteris- 
ing radiation reflected from disc, including the iron line, and 
the primary hard X-ray continuum emission. Namely, weak 
variability of the reflected component, uncorrelated with the 
variability of the primary emission, has been reported in a 
number of sources (e.g., Vaughan & Fabian 2004; see re- 
view in Fabian & Miniutti 2005). This is contrary to ex- 
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pectations from a simple geometrical model of a hard X-ray 
source located close to the reflecting disc, where a strict 
correlation between variations of the primary and reflected 
emission should be observed. 

Fabian & Vaughan (2003) first argued that such a re- 
duced variability may be explained by relativistic effects, in 
particular by light bending and focusing the primary emis- 
sion towards the accretion disc. Qualitatively, variations of 
the reflected emission should be much weaker as changes of 
the height of the X-ray source cause variations of its ob- 
served luminosity at infinity, while changes of the flux re- 
ceived by the disc (enhanced by the gravitational focusing; 
e.g., Matt et al. 1992; Martocchia & Matt 1996; Petrucci & 
Henri 1997) are much weaker. However, for a static primary 
source located on the symmetry axis, the illuminating ra- 
diation is focused into the innermost part of the disc and 
the reflected emission is subject to similar light bending as 
the primary. As a result, similar variability characterises the 
primary and reflected emission, at least for observers with 
low inclination angles. 

Miniutti et al. (2003) and Miniutti & Fabian (2004) 
have developed further this model to include rotation of the 
primary source around the axis and the resulting beaming 
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of its emission toward outer regions of the disc. Then, vari- 
ations of the reflected emission are reduced because of the 
more extended region it is produced. Predictions of their 
model are found to be consistent with observations of black- 
hole systems, e.g. by Miniutti et al. (2003), Miniutti, Fabian 
& Miller (2004), Fabian et al. (2004). However, a number of 
important effects were not systematically studied. In par- 
ticular, a specific pattern of motion of a primary source 
is assumed (including both location relative to the symme- 
try axis and azimuthal motion) but it is not discussed how 
strongly the resulting properties depend on these assump- 
tions. For example, corotation of the primary source with 
the disc is assumed by Miniutti & Fabian (2004), which is 
likely close to the disc surface, but at most approximate 
at high latitudes. As the azimuthal motion of the primary 
source relative to the disc may affect significantly the re- 
flected emission (see e.g. Reynolds & Fabian 1997), it is not 
clear whether predictions of this model are specific to the 
underlying assumptions or generic to models involving the 
light bending. 

In this paper we systematically analyse the light bend- 
ing model, concentrating on strong gravity effects. We ne- 
glect some other effects which may contribute to the original 
variability problem, for example, ionization of the disc sur- 
face (Nayakshin & Kazanas 2002). We focus here on the 
iron Ka line; an analysis of reflected emission including the 
Compton reflected radiation will be presented in our next 
paper. We study in details a number of geometrical sce- 
narios of the source location and motion. We point out cer- 
tain inadequacies in the original computations of Miniutti 
& Fabian (2004), and how they influence their quantitative 
results. We find also a novel scenario in which reduction of 
the line variability follows directly from properties of photon 
transfer in the Kerr metric. Namely, we find that a source lo- 
cated close to the disc surface, with a varying radial distance 
from a Kerr black hole, gives rise to both an approximately 
constant illumination of the surrounding disc and a very 
strong variability of the primary emission observed at infin- 
ity. Some aspects of variability in similar models have been 
considered recently e.g. by Czerny et al. (2004) and Pechacek 
et al. (2005), however, these studies considered only emis- 
sion emerging locally from the region under the source. On 
the other hand, we find that transfer of primary emission 
to more distant regions of the disc is crucial for variability 
effects in such a model. 

We concentrate on low inclination objects, with appli- 
cation to Seyfert 1 galaxies in mind. Obviously, the Doppler 
shifts are more pronounced at high inclinations, but obser- 
vational studies of high inclination objects are less advanced, 
either because they are obscured (as Seyfert 2 galaxies), or 
because the studies on dynamical time scale are not possi- 
ble (as in stellar mass black hole systems). Furthermore, we 
consider only emission averaged over at least an orbital pe- 
riod. Effects resulting from varying azimuthal location of an 
off-axis source, with respect to observer, are studied, e.g., in 
Ruszkowski (2000), Yu & Lu (2000) and Goyder & Lasenby 
(2004). 

In Section 3 we analyse various physical effects relevant 
to formation of relativistic line profiles and variability ef- 
fects; in Section 4 we apply our results to a Seyfert 1 galaxy 
MCG-6-30-15 which is the best studied object with clear 
signatures of strong gravity effects. 



2 MODEL DESCRIPTION 

We consider an accretion disc, surrounding a Kerr black 
hole, irradiated by X-rays emitted from an isotropic point 
source (hereafter referred to as the source of primary emis- 
sion). The black hole is characterised by its mass, M, and 
angular momentum, J. We use the Boyer-Lindquist (BL) 
coordinate system x 1 = (t, R,8, (f>). We also make use of lo- 
cally non-rotating (LNR) frames (Bardeen, Press & Teukol- 
sky 1972), which are particularly convenient for studying 
processes in the Kerr geometry (in the Schwarzschild met- 
ric the LNR frame is equivalent to the rest frame of a local 
static observer). 

The following dimensionless parameters are used in the 
paper 

R J A<f> 

r = -^, a — „ , , , Q——, (1) 
R e ' cR g M di 

where R g — GM/c 2 is the gravitational radius and t — 
ct/Rg. In most cases we compute spectra for a = 0.998. To 
illustrate effects which are uniquely related to black hole 
rotation, we compare in some cases spectra for a — 0.998 
with spectra for corresponding models in the Schwarzschild 
metric (a = 0). 

We assume that a geometrically thin, neutral, optically 
thick disc is located in the equatorial plane of the Kerr ge- 
ometry. The radial coordinate in the disc plane is denoted 
by rd. For distances larger than the radius of the marginally 
stable circular orbit, r ms (= 1.23 and 6 for a = 0.998 and 
a = 0, respectively) we assume a circular motion of matter 
forming the disc, with Keplerian angular velocity (Bardeen 
et al. 1972), 

fi K (r d ) = l —. (2) 

a + r d > 

We take into account emission from matter extending 
down to the event horizon. In most cases (except for model 
So, see below) we assume that matter free falls within r ms . 
We strictly follow Cunningham (1975) in description of ve- 
locity field at rd < r ms in models assuming the free fall. In 
most cases we assume the outer radius of the emitting re- 
gion in the disc r out = 600, but we compute also spectra for 
smaller r out for comparison with previous studies. 

We assume that the primary X-ray emission is gen- 
erated isotropically in the source rest frame. We neglect 
transversal or radial motion of the source, u e = and 
u r = 0, but take into account azimuthal motion. The mo- 
tion is characterised by V = v /c, where is the azimuthal 
velocity with respect to the LNR frame. To illustrate influ- 
ence of V we compare models with V = 0, which allows to 
separate effects due to gravity from kinematic effects, and 
V = 0.5, which is a typical value of orbital velocity in the 
innermost regions. 

Inclination of the rotation axis of the black hole to the 
line of sight is given by /x ODS = cos o bs- Most results (except 
for these shown in Fig. I12|l presented in this paper are for 
Mobs = 0.85, which is typical for objects presumably affected 
by strong gravity (in particular, MCG-6-30-15 and Narrow 
Line Seyfert 1 galaxies). 
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2.1 Monte Carlo 

We use a Monte Carlo method, involving a fully general 
relativistic (GR) treatment of photon transfer in the Kerr 
space-time, to find spectra of primary emission and Fe line 
observed by a distant observer. A large number of photons 
are generated from the primary source with isotropic distri- 
bution of initial directions in the source rest frame. Lorentz 
transformation, with relative velocity V, yields photon en- 
ergy, E\ n , and momentum components (in the BL coordinate 
directions) , pr r ) , pre) and prp, , in the LNR frame. These, in 
turn, yield constants of motion (cf. eq. (10) in Niedzwiecki 
(2005)) 

/ AE\ 1/2 „ (A^ 2 
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where E lrL { is the photon energy at infinity, Q is the Carter's 
constant, L is the component of angular momentum parallel 
to the black hole rotation axis, and 



A 



r 2 — 2r + a 2 , 



2.2 2 , 
r + a cos ( 



A = (r 2 + a 2 ) 2 - a 2 A sin 2 9, u = 2ar/A (6) 

A photon trajectory is fully determined by A and rj. 

For each photon, equations of motion are solved to find 
whether the photon crosses the event horizon, hits the disc 
surface or escapes directly to distant observer. Photon emis- 
sion rate, yielding normalisation of directly observed and il- 
luminating fluxes (see below), is corrected by time dilation 
factor, 



9t 



2\l/2 



(7) 



where dt s is the time interval in the source rest frame. 

For a photon hitting the disc surface, at ra, we find 
incidence angle and energy, Edisc, in the disc rest frame. 
Then, for Edi sc > 7.1 keV, we generate an iron Ka photon, 
with energy 6.4 keV, emerging from the disc. The relative 
weight of the Fe photon is related to the initial energy and 
direction of an incident photon by the quasi-analytic for- 
mula [eqs. (4-6)] from George & Fabian (1991). We modify 
the original formula by multiplying it by a factor of 1.3 to 
account for different elemental abundances (we use abun- 
dances of Anders & Grevesse 1989, while George & Fabian 
assumed abundances of Anders & Ebihara 1982). We have 
verified the factor of 1.3 using the Monte Carlo code of Zycki 
& Czerny (1994). Initial direction of Fe Ka photons is gen- 
erated with uniform distribution of azimuthal angle but we 
consider various distributions of the polar emission angle, 
jtiem = cos# em . In most cases we assume a limb darkening 
in electron scattering limit, 7(/i C m) cc 1 + 2.06/i cm but in 
S ect ions 1 3 . 7 1 and l4l we take into account other angular laws. 

For photons hitting the disc we perform also simulation 
of Compton reflection. Then, we solve equations of motion 
in the Kerr metric for both the Fe Ka and Compton re- 
flected photons emerging from r^. The detailed analysis of 
the spectral component formed by the latter, reaching a dis- 
tant observer, will be presented in our next paper. Here we 



take into account effects due to reflected photons which re- 
turn to the disc giving rise to fluorescence (Section 13.21) . 
Furthermore, in Section ^. 4l we take into account the Comp- 
ton reflected component for computing the equivalent width 
of the Fe line. 



2.2 The models 

We study properties of the iron line, in particular effects 
related to azimuthal motion and returning radiation, for 
two extreme locations of the primary source, represented 
by models A and S, defined below, with the source located 
close to the symmetry axis and close to disc surface, respec- 
tively. In these models, location of the primary source is 
given by its BL coordinates, r s and 9 S . We illustrate effects 
resulting from varying distance of the source by changing 
r s at constant 9 S , which seems appropriate for the inner- 
most part of the flow (e.g. numerical simulations typically 
show conical coronae or outflows rather than plane parallel 
structures). We consider also a recently popular variability 
model with position of the primary source parametrised by 
h s = r s cos (9 S , defined as model C below, in which the polar 
position of the source covers the whole range of values of 9 a . 
Specifically, model C is equivalent to model S for small h s 
and to model A for high hs. 



2.2.1 Model A 

In model A we consider a primary source located close to 
the axis. All model A spectra in this paper correspond to 
9 S = 0.05 rad. Our model A with V = is virtually equiva- 
lent to models with on-axis source, e.g. Matt et al. (1992), 
Martocchia & Matt (1996), Petrucci & Henri (1997), Mar- 
tocchia, Karas & Matt (2000). Our assumption of location 
being slightly displaced from the axis allows us to study the 
impact of azimuthal motion in the polar region. 



2.2.2 ModelS, S K , So 

In model S we consider a source located close to the disc 
surface. All model S spectra shown in this paper correspond 
to 6* s = 1.5 rad, i.e. h s = 0.07r s . We discuss the dependence 
of model S on V, but in Section|4]we focus on a scenario with 
the primary source rigidly coupled to the underlying disc 
and corotating with Keplerian velocity, Qk('-s)- This specific 
case is denoted as model Sk- For a = 0.998 we consider only 
r s > r ms . All variability effects in this model are studied for 
varying distance, r s , and a constant polar angle # s . 

We also introduce another model, So, for a non-rotating 
black hole, a = 0, but assuming a quasi-Keplerian velocity 
of accreting matter in the plunging region below r ms : v r — 
and = 0.5c (instead of the usual free fall assumption). 
Such a velocity field is very similar to the velocity field of 
an accretion disc around an a — 0.998 black hole, in the 
region 1.25 < r < 6 (where vk{t = 3.5) = 0.5c; vk increases 
from 0.4c at r = 6 to 0.57c at r < 2). This model can 
thus serve as a comparison model to study effects specific 
to space-time around rapidly rotating black hole, which are 
significant for r s < 6. We emphasise that we do not consider 
this as a likely scenario, we are not aware of any results sup- 
porting such a rotating flow within r ms [although e.g. Agol 
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& Krolik (2002) argue that distant regions of the flow may 
be linked through magnetic tension, which presumably could 
affect motion of the flow in the plunging region]. We consider 
model So only to separate effects due to gravitational energy 
shift and bending of trajectories from kinematic effects. 



2.2.3 Model C 

Following Miniutti & Fabian (2004) we consider a model 
with a cylindrical-like motion of the primary source, referred 
to as model C, with the source located at constant projected 
radial distance p s = r s sin# s and changing the height, h s , 
above the disc surface. Furthermore, model C assumes that 
at each h a the source has the same angular velocity, equal 
to the Keplerian velocity at p s , Hk(Ps). This yields the fol- 
lowing velocity of the source rest frame relative to the LNR 
frame (Bardeen et al. 1972) 

All spectra for model C shown in this paper are for p s — 2. 

Note that equation |8| yields, for constant angular 
velocity, a peculiar dependence of V on source location. 
Namely, V increases with increasing h s at small heights, 
achieves a maximum value at h s ~ p s and then decreases. 
E.g. for p s = 2: V(h s = 0) = 0.57, V(h s = 2) = 0.69 and 
V(h s = 20) = 0.55. These changes of V, with h B , appear to 
be crucial for properties of model C. To illustrate impact of 
kinematic assumptions, we show in Section \3 . 5 1 variability in 
a model with primary source changing h B at constant p s , as 
in model C, but with constant angular momentum equal to 
angular momentum on a Keplerian orbit in the equatorial 
plane at p., l K [yielding fi = (g** ~ l K g**) / (g u - l K g**)]. 

Suebsewong et al. (2006) study effects of changing p s at 
constant h s , keeping parametrisation of the source motion 
introduced by Miniutti & Fabian (2004). Again, we note 
that for small h s , which are most relevant for such scenario, 
the assumed pattern of motion yields velocity field strongly 
affecting the derived properties. E.g., for h s = 0.5, V = 0.6 — 
0.8 for p s < 2, i.e. achieves values much larger than those 
of a Keplerian motion in the equatorial plane. Furthermore, 
results of both Miniutti & Fabian (2004) and Suebsewong 
et al. (2006) for large p s or h s are strongly affected by the 
assumed value of r out — 100, see Section \3. 5. 3 1 



2.3 Intrinsic luminosity 

In all models we assume that the primary emission has a 
power-law spectrum, with a photon spectral index F = 2, 
and an exponential cut-off at E cut = 150 keV [this cut- 
off energy is (marginally) relevant for the fluorescent line 
only due to effect of returning radiation] . Units for spectra, 
presented in Section 3, are arbitrary but all have the same 
normalisation, with intrinsic luminosity, Lo, yielding the en- 
ergy flux of primary emission, per unit solid angle in the rest 
frame of the source, Fe, normalised to unity at 1 keV. As- 
sumption of constant intrinsic luminosity is released in two 
cases: (i) in model Sjp~ which adopts kinematic assumption 
of model Sk, but assumes intrinsic luminosity changing with 
r s according to dissipation rate per unit area in a Keplerian 
disc (Page & Thorne 1974); (ii) in models C, Sk and So, 



considered in Section [33] which adopt assumptions of mod- 
els C, Sk and So, respectively, but neglect time dilation at 
primary source, i.e. assume the source intrinsic luminosity 
= L /g t . 

We present spectra averaged over azimuthal angle. Al- 
though in models defined above the axial symmetry is bro- 
ken, the orbital periods are much shorter (with one ex- 
ception, see below) than the timescale, 10 ks, over which 
variability properties are investigated in the highest quality 
data. Note that, for high a, azimuthal averaging is justified 
also for a source with V — 0, which rotates with respect to a 
distant observer with angular velocity to, comparable to Qk 
at small r, e.g. u ~ 0.2f^K at r = 4.5. 

For a black hole with M = 10 7 Mq the orbital period 
exceeds 10 ks in model Sk with r s > 10. Indeed, narrow, 
slightly shifted lines observed is some objects with such large 
masses are considered to arise from localised spots on the 
disc (e.g. Dovciak et al. 2004b) in a model equivalent to our 
Sk with large r s . Note, however, that this range of parame- 
ters is not relevant to effects considered in our paper. In par- 
ticular, GR effects are weak and effects of black hole rotation 
are completely negligible for such distances. We analyse vari- 
ability effects averaged over an orbital period for such large 
r a for potential applicability to lower mass systems. 

Finally, we ignore time delays between the primary 
and reprocessed radiation, related to light-travel times be- 
tween the source and various parts of the disc. For a central 
source illuminating the ring of radius Rd(^> Rg), these delay 
times are between (1 — sin 9 b B )Rd/c and (1 + sin 6 h B )Rd/c. 
For o bs = 30°, the minimum and maximum delay is 
2.5 x lO~ 6 r d M/M s and 7.5 x 10~ 6 r d M/M Q s, respec- 
tively. Then, for M = 10 7 Af Q and r d > 150, the delays 
exceed typical times, > 10 ks, for high quality spectra. 

We emphasise that in most cases these delays are not 
important to variability effects studied here, regardless of 
the value of M, as the major part of reprocessed radiation 
originates from rd < 10 and the fraction of Fe photons com- 
ing from rd > 100 is smaller than 1 per cent. Significant 
irradiation of disc at r<j > 100 occurs only for hs > 10 in 
model C. In this case the delays may affect variability for 
large M, we discuss this issue further in Section 13.5.31 We 
choose r out = 600 because for this value properties of nei- 
ther of our models are affected by neglecting reflection from 
distances larger than r out . Even in the most extreme case 
of model C with h s = 20 contribution from rd > 600 would 
not exceed 1 per cent. However, in other models a smaller 
value, e.g. r out = 100, could be safely assumed. 



3 RESULTS 

3.1 Radial emissivity of Fe Kq line. 

Fig. [T] shows radial illumination profiles, e p h, and corre- 
sponding radial emissivities of Fe photons, £f c , resulting 
from illumination of the disc surface by a point source. The 
e p h and £F e are given in number of photons, per unit time and 
unit area of the disc, hitting the disc surface and emerging 
from the disc, respectively, at distance rd. All radial profiles 
are averaged over the azimuthal angle, <f>. 

Note that the illuminating flux as well as photon emis- 
sion rate are equal when expressed both in the BL coordinate 
frame and in the disc rest frame 
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Figure 1. Radial illumination (e p h) and Fe Kct emissivity (eFe) 
profiles, averaged over azimuthal angle, for a disc illuminated by 
a point source, with V = 0: (a) model A with r s = 6; (b) model 
A with r s = 2; (c) model S with r B = 6. The dotted (black) and 
dot-dashed (red) curves show e pll in the Kerr (with a = 0.998) and 
flat space-time, respectively. The solid (magenta) curves in panels 
(a) and (c), approximately coincident with the dotted curves at 
r'd > 2, show for the Schwarzschild space-time. The dashed 
(blue) curves show ep c for the Kerr metric (and photon spectral 
index T = 2). The illuminating photon flux is normalised by to- 
tal number of photons emitted from the primary source. The ep c 
profiles are normalised to dotted curves at large r d . Deviations 
between dot-dashed and dotted curves illustrate impact of GR ef- 
fects on the radial illumination pattern. Deviations between dot- 
ted and dashed curves are primarily due to blueshift of photons 
in the disc rest frame. 



P h 



dtdS dt d dS d 



(9) 



where 7V p h is the number of photons, dt and dtd is the 
time interval in BL and disc rest frame, respectively, and 
dS and dSd is the surface area element in the BL and 
disc rest frame, respectively. Namely, the following relations 
dS d = 7(A/A) 1/2 r- 1 d5' (e.g. Petrucci & Henri 1997) and 
dt d = 7~ 1 r(A J A4) 1/2 di (cf. equation 0; 7 is the Lorentz 
factor corresponding to Kepler velocity, vk, of disc element 
with respect to the LNR frame) imply that both GR and 
special relativistic (SR) effects cancel out. Note also that 
(Lorentz) transformation between the disc rest frame and 
the LNR frame (labelled by 'In' subscript) includes addi- 
tional time retardation terms, dSin = 7 _1 (1 -px/c)" dSd 
and dti n = (1 — /iuk/ 'c)7<iid (where p is cosine of the an- 
gle, in disc rest frame, between the direction of motion of a 
surface element and the direction of emitted/incident pho- 
ton), which again cancel out (Mathews 1982, Laor & Net- 
zer 1990), yielding dti n d,Si n = dtdd<5d (the same conclusion 
follows trivially from Lorentz invariance of the space-time 



volume between the LNR frame and the disc rest frame; 
Petrucci & Henri 1997). 

For clear illustration of effects resulting from light bend- 
ing, we compare in Fig. [T] the e p h profiles for a black hole 
accretion disc and a disc in flat space-time, with the same 
position of the primary source when expressed in BL and 
spherical coordinates, respectively. The profiles for a static 
point source in flat space-time are determined by obvious 
geometrical effects involving the distance from the primary 
source to disc element and projection of the disc area, yield- 
ing 



e P h()"d) 



[r'd + h 2 + p 2 



2r dPs cos(<f>)]- 3/2 d<f>, (10) 



where p s (= r s sin 6$) is the distance of the source from the 
symmetry axis, h s (= r s cos# s ) is the height of the source 
above the disc surface and <j> is an azimuthal angle in the 
disc plane. For p s <C hs, Figs. [Ha)(b), the profile is flatter 
in the central region (at rd < h B ) and steepens to e p h oc r7 3 
at rd 3> h s . For p s > h s , Fig. HJc), the profile has a steep 
local maximum at rd ~ p s - 

In the black-hole model, light bending enhances illumi- 
nation of the central region. For a source close to the symme- 
try axis the deflection of light gives rise to e p h significantly 
steeper, in the innermost few R s , than in the flat model. 
Note, however, that only r s <2, Fig. [ljb), leads to steep il- 
lumination profile, with 9^3, where by q we denote index 
of a power-law approximating locally the radial profiles, i.e. 
e oc r^ q . For a source located close to the disc surface, Fig. 
[TJc), deflection of photon trajectories increases illumination 
of the opposite (to the source) side of the disc, enhancing 
e p h both at rd < r s and (moderately) at rd > r s . 

Dashed curves in Fig. [1] show emissivity of Fe Kq pho- 
tons for the model with a = 0.998. Within rd < 4, the £Fe 
profile is steeper than e p h- This effect is primarily due to 
the gravitational blueshift of photons hitting the disc sur- 
face, which increases the number of photons with energies 
above the photoelectric threshold. Note, however, that very 
steep (q > 4) eFe occurs only in the part of the disc extending 
within the ergosphere [the ergosphere is the region contained 
within the surface given by r crg = 1 + (1 — a 2 cos 2 #) 1//2 (w 2 
for a = 0.998 and 6 — 7r/2)] where the gravitational blueshift 
is sufficiently high, e.g., E d i BC /E s > 4 at rd = 2 and > 10 
around rd = 1.2(« r ms ), where E s is the photon energy in 
the source rest frame. 



3.2 Bending to the equatorial plane 

Trajectories of photons emitted from a source located at 
r s > 4 are bent toward the centre, which effect operates in 
the same way regardless of the black hole spin value. As a 
result, identical profiles of e p h (at rd > 2) are achieved for 
high and low value of a, see Figs. [Ha)(c), in models with 
r s > 4. For a source located within r s < 4 from a rapidly 
rotating black hole, an effect unique for the Kerr metric 
operates, namely the emission is strongly focused along the 
equatorial plane (Cunningham 1975, Dabrowski et al. 1997). 
This property of the Kerr space-time leads to significant en- 
hancement of illumination of the disc lying in its equato- 
rial plane, particularly for a source located within the er- 
gosphere. Strength of the bending decreases with increasing 
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Figure 2. The figure illustrates dependence of e p h (panels a-c) 
and ep c (panels d-f ) on azimuthal velocity of the primary source. 
Panel (a) illustrates also bending to equatorial plane in the Kerr 
metric by comparing e p h for a = and a = 0.998. Panel (d) shows 
effect of returning radiation on ep c for V = and V = 0.5. All 
e p h and £F e , except for dotted (blue) curve in panel (a), are for 
a = 0.998. Panels (a-c) show e p h resulting from direct illumination 
(i.e. neglecting returning radiation) for a source, with V = 0.5 
(solid curves, black) and V = (dashed curves, red): (a) model S 
with r B = 2.5; (b) model S with r B = 6; (c) model A with r B = 6. 
The dotted (blue) curve in panel (a) is for model S with a = 0, 
V = and r B = 2.5. (d) Model 5 with r B = 2. The solid (black) 
and dashed (red) curves are for V = 0.5 and V = 0, respectively, 
and result from direct illumination; the dotted (blue) and dot- 
dashed (green) curves are for V = 0.5 and V = 0, respectively, 
and take into account the returning radiation, (e) Model S with 
r B = 6. The solid (black) and dashed (red) curves are for V = 0.5 
and 0, respectively, (f ) The solid (black) and dashed (red) curves 
are for model A, r B = 6, with V = 0.5 and V = 0, respectively. 
The dotted (blue) and dot-dashed (green) curves are for model 
C, h s = 8, with V ~ 0.6 [resulting from equation JH}] and V = 0, 
respectively. 

r B , but for 2 < r B < 4 this effect remains crucial, e.g. a source 
with V = located at r s = 2.5 yields e p h by a factor of 2 
higher for a — 0.998 than for a = 0, see Fig. [2ja). Crucially 
for the strength of the Fe line observed by low-inclination 
observers, bending to the equatorial plane enhances irradi- 
ation of disc at rd > r s , while bending toward the centre, 
which is a dominant effect for higher r s , enhances irradiation 
at r d < r s . 

Deflection of photon trajectories toward the equatorial 
plane affects both the primary emission and radiation re- 
flected from inner parts of the disc (see Cunningham 1976) . 
The former case is crucial for properties of model S as we 
discuss in Section 13.51 The latter results in redistribution 
of illuminating flux on the disc surface and in some cases 
yields significant increase of the Fe line flux, as we show in 
this section. 

Enhancement of Fe Kq line by returning radiation was 
studied by Dabrowski et al. (1997) in a model with a source 
corotating with the disc, consistent with our model Sk - They 
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Figure 3. Effect of returning radiation on profiles of Fe Ka line 
for /i D bs = 0.85. The dashed (blue) profiles result from direct 
illumination only, the solid (black) profiles take into account ra- 
diation returning, (a) Model S with V = 0.5 and r B = 2. (b) 
Model S with V = and r B = 2. (c) Model A with V = and 
r B = 2. 

find that the returning continuum photons increase the flux 
of the line by no more than 20 per cent. In generic agree- 
ment with Dabrowski et al. (1997), we find in model 5* with 
V — 0.5 (and similarly in model Sk), with small r B , that re- 
turning continuum photons increase ep a at r<j > r B by > 10 
per cent, see Fig. Hid), yielding enhancement of the line by 
the same order of magnitude, see Fig. (3ja) . This property 
is obvious, as (i) the reflecting albedo around the Fe Ka 
edge is about 10 per cent and (ii) transfer of photons from 
both the corotating source and the disc under the source, to 
other regions of the disc, is similar. The overall effect of the 
returning radiation has a minor impact on model involving 
primary source corotating with the disc (see also Fig. [T^b) 
below). 

However, we note that a significantly stronger enhance- 
ment of the line by returning photons may occur if the source 
does not corotate with the disc. Equations @ and @ yield 



which implies that the gravitational shift of photon energy 
between the source, r B , and the disc, ra, 

E ln (r d ) = 1 - Xuj(r d ) 

E ln {r B ) 1-Aw(r.)" V ' 

increases with decreasing A for r B > rd- Then, the average 
gravitational blueshift is higher for lower V (implying on av- 
erage lower A of emitted photons). Furthermore, for model 
S with V = 0, Doppler blueshift of illuminating photons is 
high due to velocity gradient between the source and the 
disc. Combination of the above two effects results in higher 
blueshift, for smaller V , of photons illuminating disc within 
the ergosphere and thus in the increase of magnitude of ra- 
diation reflected from < 2 - most of which returns to the 
disc, enhancing £f more efficiently than for V = 0.5. E.g. 
in model 5* with r B = 2 and V = 0, the resulting enhance- 
ment of ep c [Fig.[2ld)] and of the line flux [Fig. EJb)] exceeds 
50 per cent. Even more significant effect occurs for a source 
with a very centrally concentrated (and weak beyond a few 
R g ) direct illumination of the disc, specifically, for a source 
with V — at a small r B in model A, see Figs.|4ja) and[3lc). 
E.g. for r B = 1.6, due to contribution of returning radiation, 
the line flux increases by a factor of ~ 3. 

Enhancement of the Fe line by returning radiation is 
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Figure 4. Radial ep c profiles resulting from direct illumination 
in our basic models, (a) Model A with V = and r B = 1.6 (solid 
curve, black), 3 (dashed curve, red) and 20 (dotted curve, blue), 
(b) Model C with h B = 2 (solid curve, black), 8 (dashed curve, 
red) and 20 (dotted curve, blue), (c) Model 5k with r B = 1.8 
(solid curve, black), 2.2 (dashed curve, red), 3 (dot-dashed curve, 
green) and 4 (dotted curve, blue). The dot-dashed (green) curve 
in panel (a) takes into account contribution of returning radiation 
for r s = 1.6. 

efficient only for small r s . For r s > 3, the flux of return- 
ing radiation is in general much smaller than the directly 
illuminating flux. Note also that all effects discussed in this 
section are unique for high a models. Specifically, we find 
that for a — returning radiation enhances the Fe line flux 
at most by a few per cent. 
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Figure 5. Impact of the primary source azimuthal velocity on the 
Fe Ka line profiles for fi \> B = 0.85. Panels (a) and (c) compare 
profiles for V = 0.5 (solid, black) and V = (dashed, blue) in 
model A with r B = 6 (a) and in model S with r s = 2 (c). Panels 
(b) and (d) are for model C with h s = 20 (b) and h B = 0.5 (d); 
the solid (black) profiles are for V resulting from equation J5}, i.e. 
V = 0.55 (b) and V = 0.59 (d); the dashed (blue) profiles are for 
the same locations of primary source but with V = 0. All profiles 
include contribution of returning radiation. 



3.3 Azimuthal motion 

We point out that azimuthal motion of the primary source 
affects strongly the line flux through combination of GR 
and SR effects. The relevant SR effects are discussed by 
Reynolds & Fabian (1997). Doppler shifts and aberration 
of the primary emission influence both the number and av- 
erage direction of primary photons above the iron photo- 
electric threshold in the disc rest frame. On the other hand, 
azimuthal motion reduces slightly both the illuminating and 
observed flux due to Lorentz time dilation (by I/7 ~ 0.87 
for V = 0.5). 

Additional effects occur for a black hole accretion disc, 
most pronouncedly for small r B (< 4) and high a. With in- 
creasing V, collimation toward direction of motion (+<j) di- 
rection) increases fraction of photons bent toward the equa- 
torial plane and thus enhances e p h. For 2 < r B < 4, this effect 
is moderate, e.g. in model 5* with r B = 2.5 increase of veloc- 
ity from V = to V = 0.5 results in increase of e p h at rd > 3 
by 20 per cent, see Fig. E^a). Much more efficient increase 
of the fraction of primary photons bent to the disc - with 
increasing velocity - occurs for r B < 2. Furthermore, the SR 
collimation increases the average p^j and therefore the av- 
erage energy of primary photons increases, see equation ©, 
yielding additional enhancement of the ionizing flux of pho- 
tons hitting the disc. Combination of these effects results in 
a strong enhancement of £Fe by increasing V. E.g., in model 
S with r B — 2, increase of V from V = to V = 0.5 results 
in the increase of ep c by a factor of 2.5, Fig. Eld). 

Note that the enhanced irradiation of the disc by a 
source with higher V, discussed above, results only from the 



direct emission from the source. On the other hand, models 
with smaller V are characterised by stronger contribution of 
returning radiation, as discussed in Section 3.2. Figs.[5jc)(d) 
show example profiles of lines for V — and V ~ 0.5, in 
models with r B ~ 2. 

For r B > 4 changes of the line caused by change of V are 
more subtle, see Figs. [5ja)(b), and result mostly from the 
Doppler blueshift. As shown in Figs. HJb)(c), increase of V 
approximately does not change e p h at rd > r B (assuming that 
the photon emissivity rate in the source rest frame does not 
change). Although the SR collimation of primary emission 
into the direction of motion increases the fraction of photons 
illuminating the disc at rd > r B , this increase is balanced 
by the SR time dilation. On the other hand, increase of V 
reduces e p h in the central region of the disc, more noticeably 
for small h B (compare Fig. [2jb) and Fig. HJc)), which effect 
is, however, not important for face-on observers. Then, for 
high r s , dependence of ep c (and hence of the emitted line) 
on V is mostly due to Doppler shifts. These energy shifts are 
slightly higher for increasing h B due to higher, on average, 
velocity gradient between the source and the disc. 

In general, the above GR and SR effects enhance the il- 
luminating flux at the expense of the primary flux observed 
at low inclinations, which increases sensitivity of variabil- 
ity models (Section 3.5) on assumptions regarding the az- 
imuthal motion. As in Section 3.2, we emphasise that the 
GR effects discussed in this section are specific for models 
involving a rapidly rotating black hole. 
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Figure 6. Impact of the space-time metric on the Fe Kct line 
profiles for fj, ohs = 0.85. (a) Model A, with V = and r s = 4. The 
solid (black) and dot-dashed (green) curve is for a = 0.998 and 
a = 0, respectively. Panels (b— f) are for models S with V = 0.5: 
(b) r s = 3, a = 0.998; the dashed (red) curve shows contribution 
of photons emitted from a narrow ring (r^ = 3.5 — 4.5) under 
the source, the dotted (blue) curve shows contribution of photons 
emitted from > 7; (c) r B = 4, a = 0.998; the dashed (red) 
and dotted (blue) curve shows contribution of photons emitted 
from r\j = 4.5 — 5.5 and r<j > 7, respectively; (d) r s = 4, a = 
(free fall within r ms ); (e) r s = 5, a = (free fall within r ms ); (f) 
Model So with r s = 4. Profiles for a = in panels (a,d,e,f) are 
decomposed into contribution from the region of the disc within 
r ms (dashed, red) and beyond r ms (dotted, blue). 



3.4 The red wing 

In this section we focus on the shape of the red wing (i.e. < 6 
keV) of the line, which is formed mostly by photons emitted 
within a few ( < 6-7) i? g . We consider effects of space-time 
metric on the red wing shape and we indicate properties of 
model S which uniquely allow to investigate effects related 
to value of a. Applicability of various models to modelling 
the elongated wing observed in MCG-6-30-15 is discussed 
in Section \4. II 

Profiles emerging due to combined emission from var- 
ious ra have qualitatively similar shapes for low and high 
a, see Fig. [6ja)]. The smooth low-energy tails of such line 
profiles indicate that photons forming them come from a re- 
gion at a few R g but they do not possess features allowing 
to constrain more detailed properties of the emitting region, 
specifically the space-time metric. In particular, change of 
the angular law, J(pi om ), has at least the same influence on 
the line shape as a change of the value of a, compare Figs. 
[HI a) and 1161 Then, it has been argued (e.g. Dovciak, Karas 
& Yaqoob 2004a) that these time-averaged line profiles do 
not contain sufficient information to yield unambiguous in- 
formation about the black hole spin. Such a featureless wing 
is always produced as a result of illumination by a source 
close to the axis, i.e. in our model A and similarly in model C 
(except for very small h s , see discussion of model S below). 

Detection of photons with energies E a h s < 4 keV, i.e. 
redshifted by E ohB /E diBC < 0.6, is the key property which 
could be used to infer a high value of a (for small incli- 
nations). Presence of such photons in the Fe Ka line may 



distinguish high from low values of a, under the assump- 
tion that in the latter case the Fe photons emerge only from 
> r ms . If contribution from within r ms is taken into account, 
profiles in models with a — and a = 0.998 extend down 
to similar energies, see Fig. [6[a) . Models of the plunging re- 
gion (e.g. Young, Ross & Fabian 1998) assess full ionization 
at Vd < 5 for a = 0, ruling out contribution from that re- 
gion. However, they make several simplifying assumptions, 
most importantly, they base on properties of a test particle 
motion. 

Regardless of these details, escape probability from the 
region producing photons observed at < 4 keV is small. 
Therefore, even for a source at a small height on the symme- 
try axis, the resulting Fe flux below 4 keV is rather weak, see 
the solid profile in Fig.JfjJa), and it may be difficult to decon- 
volve this emission from continuum components. Generation 
of a pronounced wing at such energies requires a source lo- 
cated close to the disc surface in the region emitting pho- 
tons with the most extreme redshifts (see also Section [4.1[) . 
We point out that effects of space-time metric can be stud- 
ied in line profiles emerging in such cases, specifically, for a 
source located at (i) r s < 4 and (ii) close to the disc surface 
(h s < 0.1r s ). In such a case, most of fluorescent photons 
are emitted from a relatively narrow ring, with the width 
Ar<j < 1, under the primary source, see Fig. [4jc) . Emission 
from this ring has a well defined gravitational redshift as well 
as Doppler shifts. For /i Q b s = 0.85, magnitude of the grav- 
itational redshift (-Emf/^LN < 0.7 for rd < 4) exceeds the 
Doppler blueshift and photons emitted from rd ~ r a form 
a clear feature, indicated by dashed curves in Figs. HHb,c), 
with a pronounced blue horn (at 4-5.5 keV, depending on 
r s ) and a strongly suppressed red horn. Note also that the 
overall profiles always have a strong blue peak around 6.5 
keV, resulting from emission from larger ( > 7R g ) distances. 
While energy of the latter allows to constrain the inclination 
angle, the energies and heights of the horns from rd ~ r s al- 
low to study effects of space-time metric. Namely, equation 
([21) implies that the gravitational redshift of photons emit- 
ted from rd ^ 4 depends on a noticeably; e.g. for /i t, a = 0.85 
and rd = 4, photons are redshifted by i?i n f/i?LN = 0.71 for 
a = 0, while for a = 0.998 their redshifts are in the range 
E inf /E LN = 0.66 - 0.77 (for a = 0.998 the redshift depends 
on photon initial direction via pu))- Then, ta 10 per cent 
difference in energies of the horns (formed at rd ~ r s ) may 
be expected between low and high a. Furthermore, for pho- 
tons originating from rd < 4, the space-time metric affects 
trajectories, which is reflected in heights of the horns. 

Indeed, such differences between model S with a — 
0.998 and model So, both with r s = 4, are clear in Fig. 
[|)Jc,f). Note that these differences result only from GR ef- 
fects - in both cases the disc under the source has vk — 0.5c. 
In Fig. |6jd,e) we show spectra for model S with a — (and 
the usual assumption of free fall within r ms ). Due to addi- 
tional Doppler redshift (v r < 0) the redshifted horns have 
lower energies and are suppressed with respect to the blue 
peak more significantly than in models So with the same r s . 

3.5 The GR models of spectral variability 

Uncorrelated variability between direct continuum and re- 
flection components (including the iron line) revealed in sev- 
eral objects is difficult to understand if we observe the same 
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Figure 8. Changes of Fe Ko line profile, for fi b s = 0.85, in 
models A and C with a = 0.998. (a) Model A, V = 0, with 
r B = 1.6 (solid, black), 2 (dashed, red), 3 (dot-dashed, green) and 
6 (dotted, blue), (b) and (c) are for model C and C, respectively, 
with h B = 2 (solid, black), 4 (dashed, red), 8 (dot-dashed, green) 
and 20 (dotted, blue). Profiles in panels (b) and (c) are rescaled 
by a factor 0.4. 
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Figure 7. Total flux in the Fe Ka line vs the primary emis- 
sion flux observed at /i D b s = 0.85. All models in panels (a) and 

(b) assume a constant intrinsic luminosity of the primary source, 
changes of the observed fluxes result from varying position of the 
source (in general increasing distance corresponds to increase of 
the fluxes), (a) The solid (black) curve is for model Sk with 
a = 0.998 and the radial distance between r B = 1.5 and r B = 100; 
the (blue) diamond is for r s = 4. The dashed (red) curve is for 
model A with a = 0.998, V = and r B = 1.6-20. The dotted 
(blue) curve is for model C with h B = 0.07-20; the (red) dia- 
monds are for h B = 1 (left) and h B = 8 (right). The heavier 
dot-dashed (green) curve is for model So and the thinner dot- 
dashed (magenta) curve is for model 5k with a = (and a free 
fall inside r ms ); both models with a = are for r B = 3-100 and 
assume V = 0.5 at r B < 6. The heavy solid (magenta) curve is for 
constant p B = 3, constant (Keplerian) angular momentum and 
h B = 0.07—20. All models in panel (a) take into account the re- 
turning radiation, (b) The solid (black) and dashed (red) curve is 
for model Sk and A, respectively, same as in panel (a). The dot- 
ted and dot-dashed curves are for the same parameters but with 
neglected effect of returning radiation, (c) The solid (black) and 
dashed (red) curve is for model S K T (i.e. for intrinsic luminos- 
ity changing with r B according to radial emissivity of a Keplerian 
disc) with a = 0.998 and a = 0, respectively. The primary source 
changes distance between r B = 1.5 (for a = 0.998) or 6 (for a = 0) 
and r B = 100 (toward the bottom). Normalisations of curves in 

(c) are chosen to yield the same maximum values of the direct 
continuum flux. 



Figure 9. Changes of Fe Ko line profile, for /i b s = 0.85, in model 
S. Panels (a) and (b) are for model 5k and Sk (neglecting the 
time dilation), respectively, with a = 0.998 and r B = 2.2 (solid, 
black), 2.5 (dashed, red), 3 (dot-dashed, green) and 4 (dotted, 
blue). Panel (c) is for model So with a = and r B = 3 (solid, 
black), 4 (dashed, red), 4.5 (dot-dashed, green) and 5 (dotted, 
blue). 



primary emission that illuminates the disc. As proposed by 
Miniutti et al. (2003) and Miniutti & Fabian (2004), vari- 
ability of both the observed continuum and line may be in- 
duced by light bending if the primary source is concentrated 
near the axis of a Kerr black hole. Moreover, they find that 
variations of both components are not correlated for some 
range of parameters. 

Following Miniutti & Fabian (2004), we compute the 
observed flux of the primary source vs. that of the Fe line, 
as the source position is changed, for a source with constant 
intrinsic luminosity, see Fig. [7|a). Changes of ep c , related 
to varying position of the source, in our basic models are 
shown in Fig. [4] the corresponding changes of line profiles, 
for model A, C and Sk, are shown in Figs. [8fa). [E|b) an d 
[9ja), respectively. 
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3.5.1 Model A 

In model A, the primary continuum flux changes by two 
orders of magnitude when the distance of the source varies 
within r s = 1.6 — 20. However, the line flux varies similarly 
in this model. At small r s , directly observed primary flux is 
significantly reduced while the illuminating flux enhanced by 
the light bending. The enhanced illumination of the disc is, 
however, concentrated in the innermost region (at rd</i s , 
see Fig. [TJ and the reflected emission is subject to similar 
bending (due to which most of this emission is lost into black 
hole) as the primary emission. As a result, the line observed 
at low inclinations [Fig. [H£a)] roughly follows changes of the 
continuum. Variations of the line flux are reduced by a factor 
of ~ 2 with respect to the continuum flux, primarily due to 
returning radiation, see Fig. 0b) . 

3.5.2 Models S 

The most significant reduction of variability of the line flux 
seems to characterise model Sk, with a = 0.998 and r s < 4, 
where the line flux changes by a factor of 7 while the primary 
continuum flux changes by two orders of magnitude. On the 
other hand, in model So (with the same kinematic properties 
but a = 0) both components are linearly correlated, while in 
model Sk with a — (involving free fall within r ms ) inward 
acceleration in the plunging region results in increasing de- 
cline of the line flux; see the dot-dashed curves in Fig. [7J a )- 
For r s > 6 model 5k yields the same relation between the 
line and continuum fluxes (regardless of a) as model A with 
V = 0. 

We stress that the reduced variability of the line flux, 
for small r s and high a, is directly related to properties of 
the Kerr metric. Namely, decreasing r s results in stronger 
bending to equatorial plane and this in turn gives rise to 
enhanced illumination of the disc surface at > 6, where 
photons forming the blue peak around 6.5 keV come from, 
see Fig. H){b,c). As a result the line profiles, Fig. HJa), may 
be decomposed into two components: (i) the redshifted emis- 
sion from the region under the source, which changes sim- 
ilarly to the primary emission; and (ii) the blue peak with 
much lower variability. 

The distinct properties of the Kerr metric, underlying 
reduction of the line variability in model S, are clearly il- 
lustrated in models neglecting the effect of time dilation at 
the primary source, i.e. with the source luminosity = Lo/gt, 
where g t is given by equation (0; note that the time dila- 
tion affects in the same way both the directly observed and 
illuminating flux, thus increasing magnitude of variability 
of both the primary and reflected components, but does not 
change the relative variability of these two components. 

The eFe(r) profiles in models neglecting time dilation 
are shown in Fig. 1101 In model Sk with a = 0.998, focusing 
of primary emission toward the disc yields constant £Fe(f ) in 
the region where photons forming the blue peak originate. 
This, in turn, results in identical line profiles at E > 6 keV 
(Fig. Ob)), while the primary emission flux changes (simi- 
larly to the line flux at < 6 keV). A similar neglect of time 
dilation in model So does not result in any range of ra with 
similarly constant eFe, Fig. HOf b); in this case flux in the 
blue peak is strictly correlated with the primary flux, see 
Fig. He). 




Figure 10. Changes of eF c (r) resulting from the change of dis- 
tance of primary source in models S (with neglected time dila- 
tion). Panels (a) and (c) show impact of azimuthal velocity pat- 
tern, Keplerian vs constant with V = 0.5, for a = 0.998. Panels 
(b) and (c) show effect of black hole spin, a = vs 0.998, for con- 
stant V = 0.5. (a) The solid (black), dashed (red), dot-dashed 
(green) and dotted (blue) curves are for model Sk with a = 0.998 
and r s = 1.8, 2.2, 3 and 4. (b) The solid (black), dashed (red) and 
dotted (blue) curves are for model Sq with V = 0.5, and r s = 3, 
4 and 5. (c) The solid (black), dashed (red) and dotted (blue) 
curves are for model S with a = 0.998, V = 0.5 and r s = 1.6, 2 
and 3. 



Note also that V increases with decreasing r s in model 
Sk, which seems to be the most likely for a source close to 
the disc surface, as it may be expected to be rigidly coupled 
to the accretion disc by magnetic fields. In general, the ob- 
served line and continuum fluxes are sensitive to assumed 
V, see Section 3.3 and discussion of model C below. We 
emphasise, however, that the increase of V in model Sk is 
not crucial for reduction of variability of the line. For high 
a and small r s , bending to the disc plane is a robust effect, 
exceeding in magnitude other relevant effects. In particular, 
approximate constancy of eFe(r) around = 10 is achieved 
also for fixed V, see Fig. HOf c). 

Our simplified description of generation of X-rays by a 
point source with constant intrinsic luminosity allows to il- 
lustrate various GR and SR effects, as in panels (a) and (b) 
of Fig. [7] Combination of these effects with more realistic 
modelling of the X-ray source may result in variety of ad- 
ditional properties. A particularly interesting case of model 
S K T , with intrinsic luminosity proportional to the energy 
release rate in a Keplerian disc, is shown in Fig. [7J C )- Ef- 
fects due to bending to the equatorial plane are again clear 
in comparison of models involving a = and a — 0.998. In 
the former case, both the line and continuum fluxes increase 
with decreasing r s (> 10) due to increase of the dissipation 
rate, however for r s < 10 the fluxes decrease, due to decrease 
of both the dissipation rate and photon escape probability, 
following approximately the same relation between the con- 
tinuum and line flux as at higher r s . On the other hand, for 
a = 0.998 the primary emission flux achieves the maximum 
value around r s = 3 and than decreases by a factor of 3 
with the distance decreasing to r s = 1.8, while the Fe flux 
remains almost constant. 
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Figure 11. The Fc Ka line profiles, for /x t, s = 0.85, in model C 
with h B = 4 (a), 8 (b) and 20 (c). The solid (black) profiles are 
for r ou t = 600 and the dashed (red) are for r ou t = 100. 



3.5.3 Model C 

Miniutti et al. (2003) and Miniutti & Fabian (2004) show 
that with modifications of the model with source close to the 
axis (see definition of our model C), variability of the line 
flux is strongly reduced over some range of heights. Specifi- 
cally, they identify three different regimes in which the iron 
line is correlated (at small h s ), almost constant (at interme- 
diate h s ) and anti-correlated (at high h s ) with respect to the 
direct continuum. For fi h s = 0.85 and p s = 2, these regimes 
correspond approximately to h s = 1 — 4, h s — 4 — 13 and 
h s > 13, respectively. 

The dotted curve in Fig. [7|a) shows the iron line flux 
[changes of the line profile are shown in Fig. [SJb)] as a func- 
tion of primary continuum flux for model C. In agreement 
with Miniutti & Fabian (2004) we find that the primary 
continuum changes by a factor of 20 for the height varying 
between h s = 1 and 20. On the other hand, we find differ- 
ent dependence of the line flux on h s , in particular we do 
not note a decrease of the line flux at large heights. We sus- 
pect that this discrepancy results primarily from different 
values of the outer radius of the disc, r ou t = 100 assumed by 
Miniutti & Fabian (2004) and r out = 600 in this paper. An 
additional, small difference involves angular law of Fe emis- 
sion (assumed to be locally isotropic in Miniutti & Fabian 
2004), which yields a moderate (<20 per cent, see Fig. I17[) 
change of quantitative results. On the other hand, value of 
r ut appears to be crucial for qualitative predictions of the 
model. 

The finite extension of the outer radius of the disc re- 
sults in the presence of a central dip in the line profile, which 
is not important for the total flux of the line for sufficiently 
large r out . However, at h s > 10, the solid angle subtended 
by a disc with r out = 100 is significantly smaller than, e.g. 
for r out — 1000 (see Matt et al. 1992). Moreover, illumina- 
tion is shifted to outer regions of the disc in model C due 
to azimuthal motion. As a result, the line flux is increas- 
ingly underestimated with increasing h s , e.g. by over 30 per 
cent for h s = 20, in the model with r out = 100, see Fig. [TT] 
We conclude that for face-on observer the anti-correlation is 
an artificial effect which does not occur if sufficiently large 
r out is taken into account. On the other hand, we find a real 
anti-correlation for higher inclinations, see below. 

Obviously, illumination of distant regions, and hence 
the choice of specific value of r ou t, becomes more important 



for increasing h s . Our aim is to study effects related to GR 
or SR rather than simple geometric effects related to the 
change of the solid angle. Therefore, we assume r out = 600, 
for which reflection from ra > r out is negligible (for h s < 20), 
in our computations. Clearly, such a value is necessary for 
applications of the model to stellar-mass black hole systems 
(cf. Miniutti, Fabian & Miller 2004). For large black hole 
masses, M > 10 7 AfQ, part of reflection from distant re- 
gions of the disc would be subject to significant delays, ow- 
ing to light travel-time effects, see Section 2. In such cases, 
variations of the reprocessed radiation should be somewhat 
reduced due to mixing of radiation responding to primary 
emission from various positions of the source. On the other 
hand, for large h s , the Fe flux is close to constant even with- 
out the mixing effect. A more interesting effect, related to 
light-travel time in model C, for high values of M, would 
involve response of the line profile to changes of the source 
position. In particular, rapid increase of h s should result in 
such changes of the line profile as those predicted by Stella 
(1990). For black holes with lower masses, M<W 7 M Q , im- 
pact of time delays on variability properties depends on the 
range of heights and the time-scale over which the variabil- 
ity is analysed. Note, e.g., that for M — 6 x 10 6 Mq (which 
is within the range of masses allowed in MCG-6-30-15 by 
optical measurements, see Section|4]), the mixing effect could 
not be neglected in applications of model C with h s > 10. 

For r out = 600, taking into account the whole range of 
heights (< 20-Rg), we find the following relation between the 
line and continuum in model C. At h s < 1 the line varies by 
a factor of 2 while the primary continuum is almost constant. 
For h s = 1-8, variations of the line are reduced by a factor 
of 2 with respect to the continuum. For h s > 8, the line 
flux is close to constant. In the last case, majority of Fe 
photons come from the region of disc where light bending 
and beaming effects are small and changes of line flux result 
from varying time dilation at the primary source. Neglecting 
the time dilation we find that the line flux is constant while 
the primary flux increases by 50 per cent with increase of 
height from h s — 8 to h s = 20. Note, however, that although 
the total line flux does not change, the line profiles change 
significantly with h s [see Fig. HJc)] - in general, flux at 5.8 
keV remains approximately constant but the wing below 
this energy weakens while the blue peak strengthens with 
the increase of h s [note the change of the distance where 
most of Fe photons come from in Fig. 21b)] . 

We emphasise that the above behaviour is not a generic 
prediction, resulting from light bending, for a model involv- 
ing a source with changing h s . It is strongly related to spe- 
cific kinematic assumptions underlying model C, see Sec- 
tion 2.2.3. Apart from impact on the line flux, discussed in 
Section 3.3, changes of V reduce variability of the primary 
continuum at small heights (where weaker light bending is 
balanced by stronger SR beaming with increasing h a ) and 
enhance its variability at h s > 2. Assuming another pattern 
of motion, we obtain different relations between the spectral 
components. E.g., for a source changing h s at constant p a , 
but with constant angular momentum equal to angular mo- 
mentum on a Keplerian orbit in the equatorial plane at p a , 
V increases monotonically with h s . In this case, the mag- 
nitude of the line flux variability is higher than that of the 
continuum at all h s , see the heavy solid curve in Fig. [TJa). 
On the other hand, setting V — in model C we find, for 
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Figure 12. Same as Fig. 0a) but for an observer inclination 
A'obs = 0.45. Solid (black) curve is for model Sk, dotted (blue) for 
model C, dashed (red) for model A. All models assume a = 0.998. 



h s > 1, the same relation between the line and primary flux 
as in model A, i.e. shown by dashed curve in Fig. [TJa). 



3.5-4 Summary of variability models 

In summary, increased illumination of the innermost disc, 
within a few R g , resulting from deflection of light toward 
the centre, is not sufficient to yield reduced variability of re- 
flection component for small inclination angles. The major 
constraint results from the fact that reflection from this in- 
nermost region is not observed at small inclinations. Then, 
reduction of observed variability requires enhanced illumi- 
nation of disc beyond a few R s , which can be achieved in 
the Kerr metric for small r s as a result of bending toward 
the equatorial plane. Interestingly, the same mechanism re- 
duces variations of the line for a source at small heights at 
the axis, through the effect of returning radiation, as this 
radiation illuminates the disc at distances larger than the 
source height. Alternatively, a more extended illumination 
may result from the SR beaming and then the reduced vari- 
ability can be achieved for appropriate pattern of motion, 
e.g. for V decreasing with increasing h s . 

On the other hand, the above constraint does not ap- 
ply to systems with higher inclinations (cos# b a < 0.6) and 
for such systems the line flux varies with much smaller am- 
plitude than the continuum flux in all models, see Fig. 1121 
Moreover, anticorrelation between the line and continuum 
fluxes observed at such 9 b B is a generic property of all mod- 
els at large r B or h s . The decrease of the line flux, with 
increasing r s , is a purely SR effect, resulting from (i) an- 
gular distributions of the line flux, J ? (/i em ), in the disc rest 
frame, which declines with decreasing ^t om ; (ii) increase of 
the distance, r^, where most Fe photons come from with in- 
creasing r s ; and (iii) decrease of the orbital velocity of disc, 
with increasing ra, and thus less efficient beaming to high 
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Figure 13. The rms spectra for Fe K« line observed at ^ D b s = 
0.85. (a) Solid (black) curve is for model 5k with the radial 
distance changing between r s = 1.5 and r s = 100. Dashed (red) 
curve is for model C with h B = 0.07-20. (b) Solid (black) and 
dashed (red) curve is for model A with V = and h B = 1.6—3 and 
h s = 1.6-20, respectively. All models assume a = 0.998. 
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Figure 14. The rms spectra for Fe Ka line (top panels) and 
superposition of Fe line and primary continuum (bottom panels) 
observed at /i b a = 0.85. (a,c) The dashed (red) and solid (black) 
curves are for model Sk , a, = 0.998, with r s = 1.5-3 and r B = 2.2- 
4, respectively. (b,d) The dotted (blue) curves are for model C 
with a = 0.998 and h B = 0.07-1. The solid (black) curves are for 
model S with a = and r s = 3-5. The dashed (red) curve in panel 
(d) is for model C with h B = 0.07-20 (i.e. the same parameters 
as for the dashed curve in Fig. I13f a~)). Some spectra in panels (c) 
and (d) are shifted vertically for clarity. 
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3.6 The fractional variability 

We point out that model S with high a is uniquely capable 
to explain suppressed variability around 6.5 keV in energy- 
dependent patterns of spectral variability. We focus here on 
effects related to changes of the Fe line profile. A more de- 
tailed analysis of spectral variability, including the Compton 
reflection component, will be presented in our next paper. 

The fractional variability of observed spectra, as a func- 
tion of energy, is quantitatively characterised by the root 
mean square (rms) variability function. In models relating 
variability to change of position of the primary source, we 
determine the rms spectra according to the following defini- 
tion 



<E) = 



<f(E)>\ 



E 



ME)- < f(E) >] 2 
TV — 1 



(13) 



where the summation is performed over TV positions of the 
primary source, fi(E) is the photon flux in the energy 
band, E, corresponding to i-th position of the source and 
< f(E) > is the average photon flux in this band for all N 
positions. 

In Figs. I13l and ll4f a.b') we show the rms spectra for the 
line and in Figs. |14f c.d) for superposition of the line and 
the primary power-law. Taking into account the full ranges 
of source positions considered in our basic models - i.e. ra- 
dial distances r s — 1.5-100 and r s — 1.6-20 in models S 
and A, respectively, and the height h s — 0.07-20 in model 
C - we find that the rms spectra for the line increase toward 
low energies below 5 keV, see Fig. 1131 which effect is related 
to varying strength and extent of the red wing of the line. 
However, the line flux at < 5 keV is typically much lower 
than the primary continuum flux, therefore these variations 
of the red wing do not produce noticeable excesses in the 
rms spectra of the line plus power-law. On the other hand, 
both model A and C are characterised by reduced variabil- 
ity of the line flux around 5.5 keV (see Fig. [Hfb,c)) which 
gives rise to declines in the rms spectra both for the line 
(Fig. I13[) and for the line + power-law (dashed curve in Fig. 
1 14f V In turn, in model Sk at r s > 10, the line is formed 
in majority by photons from ra ~ r s and the energy of the 
blue peak changes with r s , giving rise to strong increase of 
rms spectrum at E > 6 keV, see Fig. I13f a). 

In most models, even these predicting the reduced vari- 
ability of the total line flux, a decline around 6.5 keV does 
not occur in the rms spectrum. We notice such suppression 
only in model Sk (and similarly S K T ) with a — 0.998 and 
r s < 4, where bending to the disc plane yields strongly re- 
duced variability of the blue peak, see discussion in Section 
13.5.21 The detailed shape of resulting signal in the rms spec- 
trum depends on the assumed range of r s , see Fig. I14f a.c). 
The range of distances extending deeper into the ergosphere 
corresponds to stronger suppression around 6.5 keV. For 
r s = 2-4, decline above 6 keV is less pronounced, on the 
other hand, a strongly variable and relatively strong red 
wing of the line gives rise to an excess variability around 
5.5 keV, see the solid curves in Fig. H4f a.cV 

Direct relation of these effects with properties of the 
Kerr metric is again clearly illustrated by comparing models 
So and 5k with a = 0.998. The rms spectra for model So 
with a (rather unphysical) assumption of the primary source 




Figure 15. The rms spectra for the superposition of Fe line and 
primary continuum observed at fi ^ s = 0.85, in model S^ T . The 
upper and lower pair of curves is for r B = 1.7—2.2 and r s = 1.8 — 3, 
respectively. The solid (black) and dashed (red) curves are for 
I(p-em) oc 1 + 2.06/x om . The dot-dashed (green) and dotted (blue) 
curves are for I(p eln ) oc log(l + l/(i em ). 



changing the radial distance at r s < 6 are shown by the solid 
curves in Figs. I14f b.d). For a — 0, variable red wing gives 
rise to excess at 5.5 keV, however, the blue peak changes 
similarly to the primary continuum and the rms spectrum 
is flat above 6 keV. 

We did not find significant differences between predic- 
tions of models Sk and S K T , see Fig. 1151 On the other hand, 
the quantitative properties are slightly dependent on the an- 
gular emissivity law, see next section. 

Finally, note that source located close to the disc surface 
(as required to model a pronounced red wing, see Section 
[3]|, but changing the height rather than the radial distance, 
gives rise to spectral variability significantly different than 
model Sk- In particular, in model C with small h s , the line 
varies much more significantly than the primary emission 
[see dotted curve in Fig. [7fa)] yielding an extended (over 
AE m 3 keV) excess in the rms spectrum, see dotted curves 
in Fig. I14f b.d). In model A with small r s (see Figs.[8ja) and 
I13f b)) the most significant changes of the line occur close to 
the blue peak. 



3.7 Limb darkening 

As noted by Beckwith & Done (2004), for a fixed radial emis- 
sivity change of the limb darkening may strongly influence 
the line profile. Therefore, we have checked our results for 
various angular laws considered in literature. 

The limb darkening in the form 7(/i e m) oc 1 + 2.06^i C m, 
assumed in most of our calculations, approximates emission 
from a disc with ionised atmosphere. The same angular dis- 
tribution is assumed in the LAOR model (see Laor 1991), 
which is often applied in spectral analysis. Emission from 
a uniform, not ionised slab is approximated by I(fj, cnl ) oc 
log(l + 1/Vcm), Ghisellini, Haardt & Matt (1994). Emission 
from an optically thin material has I([i cln ) oc l//i cm . Fi- 
nally, a locally isotropic emission, 7(/i e m) = la, is sometimes 
considered for simplicity. 

Fig. [16] shows line profiles for various angular emission 
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Figure 16. Fe Ka line profiles observed at fi h B = 0.85 for various 
angular laws. In all panels the solid (black) profiles are for limb 
darkening in electron scattering limit, I(fJ, em ) oc 1 + 2.06/x cm . (a) 
Model C with h B = 0.4. The dotted (red) curve is for locally 
isotropic emission, J(/t em ) = 7q i t ne dashed (blue) curve is for 
I(fiem) oc log(l + 1/jUem). (b) Model C with h B = 1. The dot- 
dashed (magenta) profile is for optically thin emission, I(fi C m) oc 
1/Hem- (c) Model S with r B = 2 and V = 0.5; the dashed (blue) 
curve is for /(/lem) oc log(l + l//t C m). 
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Figure 17. Total flux in the Fe Ka line vs the primary emis- 
sion flux observed at /J. a h B = 0.85. The solid (black) and dashed 
(red) curve is for model 5k and A, respectively, with a = 0.998 
and constant intrinsic luminosity, same as in Fig.[7{a), but with 
I(p-em) oc log(l + l//i e m). The dotted (blue) and dot-dashed 
(green) curve is for model C with I(fi e m) oc 1 + 2.06/i Gm (same 
as in Fig. [7Ja)) and I(fj, C m) = respectively. 



laws. Note that change of the angular law results, in some 
cases, in a significant change of height of the blue peak, while 
the lowest energy ( < 3.5 keV) part is roughly unaffected. 
Then, the slope of the red wing - connecting the low energy 
part with the blue peak - changes, following the change of 
the height of the latter. The following effects result in such 
properties of the line. 

For /i bs = 0.85, the change of the angular law turns 
out to be unimportant for the magnitude of contribution of 
photons emitted from the disc within the ergosphere (form- 
ing the observed line at _Ei n f < 3.5 keV). In this part of the 
disc, only photons emitted into a narrow range of initial di- 
rections, with intermediate values of ^t om = 0.4 — 0.6, are ob- 
served at /i bs = 0.85. Then, flux emitted in these directions 
remains roughly unaffected by the change of angular law, 
under the condition that total energy of fluorescent photons 
emitted from unit area of the disc remains constant. 

On the other hand, the blue peak is mostly due to emis- 
sion from outer regions, where light bending is not crucial. 
Then, changes in the height of the peak simply correspond 
to change of the flux emitted to /i cm « /i b s , again under 
the condition that the total emitted energy remains con- 
stant. E.g., /(fi cm ) oc 1 + 2.06/icm yields F(n em = 0.85) 
higher by a factor of ~ 2, 1.7 and 1.2 than 7(^ om ) oc l//iiem, 
I(fiem) oc log(l + 1/jUem) and I(fiem) = Io, respectively. 
These factors are in a rough agreement with differences in 
the heights of the blue peaks in Fig. 1161 

In conclusion, we confirm that change of the angular 
law results in the change of the line profile, mostly in the 
strength of the blue peak. However, we find that these (sys- 
tematic) changes of the line do not affect the main conclu- 
sions of this paper (assuming that the same angular law 
characterises emission from the whole disc surface). Specif- 
ically, the qualitative properties of variability of the total 
flux of the line as well as the rms spectra, are negligibly af- 
fected by the choice of a particular law. Figs. I15land ll7l (in 
comparison with Fig. [7J) illustrate changes of such properties 
corresponding to the change of from limb darkened to 
limb brightened. A slightly weaker signal in the rms spec- 



trum, for the latter, appears to be the most significant effect 
resulting from such a change of I (ft). Note also that changes 
of the line profile, at E < 5 keV, resulting from the change 
of the angular law are much smaller than changes result- 
ing from significant change of h s or r B , compare Figs. [5] and 
1161 making our conclusions on modelling the extended red 
wing in the next section independent of angular emission 
law. Only assumption of J(/i cm ) oc l/fx em (Fig. fToT b)) could 
affect our conclusions more significantly, however, such an 
(optically thin) emission from an accretion disc seems very 
unlikely. 



4 MCG-6-30-15 

The profile of broad Fe line in the X-ray spectrum of Seyfert 
1 galaxy MCG-6-30-15, first clearly resolved by ASCA 
(Tanaka et al. 1995), is consistent with that predicted from 
an accretion disc inclined at 30°. In part of the ASCA ob- 
servation, when the 2-10 keV flux was 2 times lower than 
average, F2-W = 4 x 10 -11 erg/s/cm 2 , for the whole obser- 
vation, the line extended below 4 keV (Iwasawa et al. 1996; 
that low flux state was defined as the deep minimum). Pro- 
files of the Fe line, extending below 4 keV, were revealed 
also in two XMM-Newton observations (Wilms et al. 2001; 
Fabian et al. 2002). 

During the first XMM-Newton observation, in 2000, the 
average flux, -F2-10 = 2.4 x 10~ n erg/s/cm 2 , was close to 
that of the deep minimum. Then, Wilms et al. (2001) and 
Reynolds et al. (2004) argued that this observation caught 
the source in its deep minimum state. However, the source 
was highly variable, changing the flux between F2-10 = 1.5 
and 4.8 x 10 -11 erg/s/cm 2 , i.e. below the deep minimum and 
above the typical flux levels from that source, respectively. 
During the second, longer observation with XMM-Newton, 
in 2001, with the average F 2 _io = 4.1 x 10~ n erg/s/cm 2 , 
the flux again varied, by a factor of 5, between values lower 
than that of the deep minimum and higher than the typical 
average value. Using the RXTE observations Vaughan & 
Fabian (2004; see their figure 1 and related discussion) argue 
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Figure 18. In all panels the dotted (black) curve shows profile of 
the Fe line derived with parameters given in Fabian et al. (2002), 
see text, (a) The solid (blue) curve is for parameters of the best 
fit in Vaughan & Fabian (2004), i.e. q ln = 6.9, q a ut = 3, T"b r = 3.4, 
n n = 1-8. The dashed (red) curve is for a model with q = 3, rj n = 
1.2 and a = 0.998. The solid (blue) curves in panels (b) and (c) 
are for model C with h B = 8 and 4, respectively, r ou t = 100 and 
7(/i em ) = Iq (as in Miniutti et al. 2003). The dashed (red) curve in 
(b) is for model A with a = 0, V = 0, r s = 3 and r; n = 6 (= r ms ). 
The dashed (red) curve in panel (d) is for model A with r s = 1.6, 
a = 0.998 and V = 0. Contribution to this profile from the region 
of the disc within and beyond r d = 2 is shown by the heavy solid 
(magenta) and the thinner solid (blue) curve, respectively, (e) 
The solid (blue) curve is for model Sk with r B = 2. (f ) The solid 
(blue) curve is for model C with h B = 0.4. All profiles, except for 
solid in panels (b) and (c), correspond to I(fi a m) oc 1 + 2.06/iom. 

that in both cases the source was observed in its typical 
state. 

We focus below on the second XMM-Newton observa- 
tion. A phenomenological model with a broken power-law 
radial emissivity (see below) gives an excellent fit to the 
data, with \v < 0.9 (Fabian et al. 2002), therefore, we pre- 
sume that this model reproduces properties of the emitting 
region with high precision and we take parameters of the 
best fit as basic for our discussion. 

The most tightly constrained value of the black hole 
mass in MCG-6-30-15, from X-ray variability, is ~ 3 x 
10 6 Mq (in agreement with optical measurements yielding 
values between 3 and 6 x 10 6 Mq; McHardy et al. 2005). 

4.1 Time-averaged profile 

As noted in Section 13.41 detection of photons with ener- 
gies < 4 keV is considered as the evidence of rapid rotation 
of a black hole. The best studied example of a relativistic 
Fe line - with elongated wing formed by such strongly red- 
shifted photons - comes from XMM-Newton observation of 
MCG-6-30-15 (Fabian et al. 2002). This is best fitted with 
a steep (gi n = 4.8) emissivity inside r br = 6.5 and a flatter 
(g ou t = 2.5) profile beyond that, and with the inner radius of 
emitting region n n = 2; the fit was obtained with the LAOR 



model (in the X-ray data fitting code xspec), which assumes 
r ut = 400 and I(fJ, em ) oc 1 + 2.06/! cm , see Laor (1991). The 
dotted profiles in all panels of Fig. [TH] correspond to those 
best fitting parameters. 

We note that there is a degeneracy between gi n and rbr, 
in general higher gi n allows for lower rb r . Vaughan & Fabian 
(2004) fit the same XMM observation (with the data set 
extended by including the EPIC pn data) with rb r = 3.4 and 
gin = 6.9, the solid curve in Fig. llSf a). Nevertheless, both 
sets of parameters require the line flux at 4 keV to be two 
times larger than in a model assuming a single power-law 
radial emissivity with q — 3 (see Fig. I18f a); the three curves 
are normalised to yield the same flux in 5-6 keV range). In 
summary, the data seem to require a steep emissivity in the 
innermost region to explain strength of the red wing. At 
larger distances the emissivity has to flatten to produce a 
rather strong blue peak. 

Such steep, broken power-law emissivity profiles can be 
produced in models with a source positioned close to the 
symmetry axis (Fig. []Ja,b); see also Martocchia, Karas & 
Matt 2000), and indeed these models were invoked to ex- 
plain the broad line profile in MCG-6-30-15 (e.g. Miniutti et 
al. 2003) . We note that in these models the source height has 
to be rather small if the steep emissivity is to extend beyond 
the ergosphere (see also Martocchia, Matt & Karas 2002) . In 
particular, we find that in model C the steep £Fe(V) profile, 
with q > 4, is produced only in the ergosphere, r& < 2, if 
h s > 1 (Fig. |4jb)). However, photons emitted from within 
rd = 2 cannot produce the extended red wing observed in 
MCG-6-30-15, since they are redshifted by 1/g = 2-3 (for 
inclination o b s = 30°), and so they are registered at E < 3.5 
keV (their contribution is shown in Fig |18f dl). Fe lines re- 
sulting from h s > 1 do have profiles skewed towards low en- 
ergies, however they do not posses elongated red wings with 
approximately constant slope between 3.5 and 5.5 keV, as 
observed in MCG-6-30-15. For example, for h s — 3 — 8, as 
required by the variability model of Miniutti et al. (2003), 
the profile is much less pronounced than that correspond- 
ing to parameters of the best fit from Fabian et al. (2002) 
(FiglTSTb.c)): compare with fig. 5 in Miniutti et al. 2003). We 
find that to fit the profile observed in MCG-6-30-15 model 
C requires the source height h s < 1, see Fig. I18f f), 

Similar constraints are obtained in model A. Here, in 
order to produce steep emissivity (g > 4) at rd > r org , the 
source must be located at r s < 4. Explaining the red wing 
in MCG-6-30-15 requires r s < 2, with the best agreement 
achieved for r s = 1.6, Fig. I18f d). This overproduces some- 
what the flux below 3.5 keV, if ri n < 2 (hence it requires 
n n ~ 2, in agreement with Fabian et al. 2002). 

Miniutti et al. (2003) claim that their model fits the 
data for h s — 3 — 8. We note that for h s > 3 the flux below 
5 keV is actually weaker than for the model with a sin- 
gle power-law with g = 3, compare the dashed curve in Fig. 
I18f a) and the solid curve in Fig. ll8f c). Moreover, for h s — 8, 
the line profile is fully consistent with the one produced by 
emission only from rd > 6 in a model with a = 0, see Fig. 
I18f b). Then, the fits in Miniutti et al. (2003) appear to con- 
tradict results of, e.g., Fabian et al. (2002) who find that 
strong emission from rd < 6 is required. We suspect that 
this discrepancy results from modelling of the continuum 
emission. The simultaneous observation with Beppo-SAX, 
taken into account by Fabian et al. (2002), allows to con- 
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Figure 19. The figure illustrates formation of a profile for super- 
position of primary X-ray sources in model with I(ii C m) oc 
1 + 2.06/^cm- The profiles for superpositions are rescaled for clar- 
ity. The solid (blue) curve in all panels shows the superposition 
profile for r s = 1.8, 2, 2.2, 2.5 and 3 (Di). (a) The dotted (black), 
red (dashed) and dot-dashed (green) curve show individual pro- 
files for r B = 3, 2.5 and 2.2. (b) The dashed (red) curve shows the 
superposition profile for r a = 1.8, 2, 2.2, 2.4, 2.6, 2.8 and 3 (D 2 ). 
The dot-dashed (green) and dotted (black) curve show individual 
profiles for r s = 2.6 and 2.8. (c) The dot-dashed (green) curve 
shows the superposition profile for sources with r s = 1.5, 1.7, 1.8, 
2, 2.2, 2.5 and 3. The dashed (red) and dotted (black) profiles are 
for r s = 1.8 and 2, respectively. 



strain the underlying power-law continuum and the Comp- 
ton reflected component. On the other hand, Miniutti et al. 
(2003) use only data in 3-10 keV range, which poorly con- 
strains the continuum. Moreover, the fitted value of V is not 
given in their paper so it is not clear if it is consistent with 
data in broader energy range. Presumably, softer power-law 
emission may compensate the red wing. We stress, however, 
that if the model of Miniutti et al. (2003) with, e.g., h s = 8, 
could fit the data in the broader energy range, it would have 
an essential implication, namely, rotation of the black hole in 
MCG-6-30-15 would not be required by the time-averaged 
profile from the XMM-Newton observation. 

In general, enhanced illumination of the disc between 
r<j = 2 and 6, required to model the strong, elongated red 
wing, seems to require the primary emission concentrated 
close to this area of the disc, as in models A and C with small 
r s or h s . Similarly, such an enhanced illumination occurs in 
model Sk with small r a . For a single primary source, in 
model 5k, the best agreement with parameters from Fabian 
et al. (2002) is achieved for r s < 2, Fig. HHe). For larger 
r s , the model predicts profiles possessing the red tail with 
a pronounced bumpy feature (cf. Figs. [6jb,c)), which is not 
seen in the time-averaged profile of the observed line. On the 
other hand, superposition of profiles corresponding to vari- 
ous r s smooths out the shape of the red wing, see Fig. 1191 
Then, we explore below in some details a model with vari- 
ability effects resulting primarily from the change of the dis- 
tance, r a , at which dominant contribution of hard X-rays is 
produced and with the time-averaged spectrum formed by 
superposition of spectra corresponding to various r s . 

Motivated by variability effects, discussed in the next 
section, we focus on model S^ T and we constrain the range 
of distances to r s < 3, which range is relevant to modelling 
the variability. Moreover, significant contribution from r s > 3 
would yield the blue peak broader than the observed one, 
at least for limb-darkened I(fi B m), see below. The detailed 



shape of the superposition profile depends (weakly for the 
considered range of r s , see below) on relative magnitudes 
of contribution from various r s . These, in turn, result from 
both the dependence of intrinsic luminosity on r s , which is 
fixed in model Sk T , and the length of time at which emission 
is generated at various r s . The latter should correspond to 
the length of periods spent by the source in various flux 
states. 

The light curve for the 2001 observation (figure 12 from 
Vaughan & Fabian 2004) indicates that the source spends 
most of the time in intermediate flux states, with rare tran- 
sitions to the very low state and more frequent, but very 
short, periods of higher flux states. We identify the lowest 
flux states with r s < 1.8 (see Section [4.6[) and the highest 
flux states with r s ~ 3. Then, the most frequent, interme- 
diate flux states would correspond to r s ~ 2 — 2.4. These 
would yield the difference by a factor of ~ 1.5 between di- 
rect fluxes of primary emission in the intermediate and high 
flux states. 

As illustrated in Fig. 1191 the detailed shape of the line is 
weakly sensitive to specific distribution of r s , provided that 
the distribution involves a few positions within r s = 2.2 — 3 
(contributions from these larger r s are shown in panels (a) 
and (b)), so that the red wing is smoothed out. In next sec- 
tions we consider model S^ T involving r s = 1.8, 2, 2.2, 2.5 
and 3 (denoted below as Di), for which the superposition 
profile reproduces the observed line, see Fig. TOf a). Further- 
more, this distribution of r s approximates the distribution 
of flux states indicated by the light-curve, with higher fluxes 
less frequently sampled, more closely than a uniform distri- 
bution of r s (denoted as D2)- Note, however, that for the 
uniform distribution D% only a weak excess between 4 and 5 
keV occurs with respect to D\, see Fig. ll9f b). These two dis- 
tributions of r s yield also similar equivalent widths (EWtot, 
see Section 14. 4[) of the Fe line, namely 670 and 620 eV for 
Di and D2, respectively. 

Individual profiles of lines for small r s (< 2) differ 
weakly from each other as well as from the profile of the 
superposition, see Fig. I19f c). Therefore, the superposition 
profile is insensitive to details of contribution from small r s , 
e.g., the profiles for r s — 1.5 — 3 and 1.8 — 3 are almost iden- 
tical (Fig. ll9f c)). This property is more thoroughly explored 
in Section T4. 6 1 

Superposition of lines from sources located at r s < 3 
reproduces shape of the observed line at E < 5 keV (see 
Fig- I20f a.b)) for both limb-darkened and limb-brightened 
emission. On the other hand, details of modelling the line 
at higher energies depend on angular emissivity law. The 
strong blue peak, at E > 5.5 keV, has to be formed by 
photons originating mostly from the inner disc rather than 
from a distant material. Vaughan & Fabian (2004) find that 
the width of the peak requires the emitting region to be 
located at rd ~ 30. However, they use a non-relativistic 
formula for the Doppler effect, therefore that value gives 
an upper limit on the distance of that region. With a fully 
relativistic modelling, we find that the peak may be repro- 
duced with the major contribution coming from rd ~ 10. For 
7(/i em ) oc l + 2.06^i om , illumination of that region by sources 
located within r s — 1.8 — 3 is approximately sufficient to 
yield the peak with the observed strength, the solid curve 
in Fig. I20f a). On the other hand, for this I(fJ, em ), strong 
emission from r s = 4 — 10 would overproduce the line flux 
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Figure 20. The dotted (black) curves in (a) and (b) show profile 
of the Fe line with parameters of the best fit from Fabian et al. 
(2002), same as in Fig. [18] (a) The solid (blue) and dashed (red) 
curve shows the superposition profile for primary X-ray sources 
distributed in the range r B = 1.8 — 3 (-Di; same as the solid 
curves in Fig. I19[ l and r B = 1.8 — 20, respectively, in model S^ T 
with /(/icm) k 1 + 2.06/x cm . (b) The solid (blue) and dashed 
(red) curve shows superposition profile for primary X-ray sources 
distributed in the range r B = 1.8 — 3 (£>i) in model SS T with 
/(/i em ) oc log(l + 1/^tom)- The solid curve includes additional 
contribution from the source located close to the axis at r B = 
10 (model [5^ T ,.A], see text), (c) The solid and dashed curve 
shows profile for h B = 8 and 4, respectively, in model C with 
I(p-em) <x 1 + 2.06/icm- The bottom panels show the difference 
spectra corresponding to models shown in top panels, (d) The 
difference between spectra for r s = 3 and r B = 1.8 in model 
SS T with ((/ion) oc 1 + 2.06/i em . (e) Same as in (d) but for 
I(Hem) oc l°g(l + 1/Vcm). (f) The solid (black) curve shows the 
difference between spectra for h B = 8 and 4 in model C with 
/(/lem) oc l + 2.06/x em . The dotted (blue) curve is for I(ii e m) = lo 
and r out = 100 as in Miniutti et al. (2003). 



in 5-6 keV range, yielding too large width of the line core, 
the dashed curve in Fig. I20f a). 

For 7(/iem) oc log(l + l//i C m), primary sources described 
by model Sf T are not sufficient to reproduce the strength 
of the blue peak, see the dashed curve in Fig. I20f b). In 
this case, an additional primary source illuminating the disc 
mostly at r<j — 10, e.g. a source at r B ~ 10 at the symmetry 
axis, is required to reproduce the blue peak. In next sections 
we discuss some properties for the specific case of sources 



described by model S£ with r B 



- 3 and an additional 



source, described by model A with r B = 10 and V = 0, 
emitting 10 per cent of total X-ray luminosity. This model, 
denoted below as [Sk T ,A|, yields the time-averaged profile 
in very good agreement with the observed one, see the solid 
curve in Fig. [20T b'). 



4.2 Variability 

A model consisting of a strongly variable, in normalisation, 
power-law and a more constant reflection component ex- 
plains variability of MCG-6-30-15 measured by ASCA (Shih 
et al. 2002) and XMM-Newton (e.g., Fabian & Vaughan 



2003). Specifically, figure 5 in Fabian & Vaughan (2003) 
shows changes of magnitude of reflection by a factor of < 2 
corresponding to changes of the primary radiation by a fac- 
tor of 5. 

Such reduction of variability could be entirely explained 
by GR effects for a source, with constant intrinsic luminos- 
ity, described by model 5k and changing the radial distance 
within r s = 2 — 3, see the solid curve in Fig.[7Ja). In this case, 
however, variations of the line, although reduced, would be 
correlated with those of the primary radiation. On the other 
hand, Fabian & Vaughan (2003) do not reveal any such trend 
in variations of both components. Then, it seems that mod- 
elling of the observed variability should involve changes of 
the intrinsic luminosity with r B . Specifically, model S^ T with 
r B — 1.8 — 3, yielding approximately constant reprocessed 
flux for the primary flux varying by a factor of 3, appears 
to be the most feasible to explain that variability pattern. 

Additional changes, by a factor of < 2, of the intrinsic 
luminosity at a given r B would be required in model Sk T to 
account for the spread of reprocessed flux values for a spe- 
cific value of the primary flux. Alternatively, such variations 
of the reprocessed component, at fixed primary component, 
could be explained by changes of h B , as for h B < 1 such vari- 
ations of the height yield an almost constant primary flux 
and the Fe flux varying by a factor of 2, see the dotted curve 
inFig.rUa). 

The reduction effect is weakened if even a small fraction 
of X-rays is generated at larger distances, as required by 
the time-averaged profile for limb-brightened I(fi), see Fig. 
I20f b~). E.g., for model [Sk T ,.i4.], involving a constant source, 
located at r B = 10 at the axis and emitting 10 per cent of 
total X-ray luminosity, variations of the distance of a source 
described by Sk within r B = 1.8 — 3 yield the change the 
observed primary flux by a factor of 2 (at the constant Fe 
flux). 

A similar reduction, by a factor of 2, of the Fe flux 
variations is predicted in model A as well as in model C, in- 
terestingly, in both cases in the range of parameters relevant 
for modelling the time-averaged profile. Namely, in model A, 
this effect - involving enhancement by returning radiation, 
see Section r3.5.1l - occurs for small r B , for which the line with 
strong red wing is produced, cf. Fig. I18f d). In model C, re- 
duced variations of the Fe flux - basing on specific V(h B ) - 
occur for the intermediate heights, for which the model fits 
the observed profile according to Miniutti et al. (2003). On 
the other hand, neither of these models is able to explain 
some finer details of spectral variability revealed in the rms 
spectra, see below. Note that effects related to light-travel 
time are (at most) marginally important in these applica- 
tions of models A and C to MCG-6-30-15. In particular, for 
h B < 8 in model C, fraction of Fe photons from r& > 150 
is smaller than 5 per cent and the longest delay time (cor- 
responding to a small fraction of photons coming from the 
farthest, to observer, side of the disc) for photons reflected 
from r d = 150 is smaller than 3.5 ks (for M = 3 x 1O 6 M ). 

Note also that the above is not the unique variabil- 
ity pattern observed in MCG-6-30-15. Namely, Reynolds et 
al. (2004) find that during the 2000 XMM observation, for 
the range of the 2-10 keV flux between F2-10 = 2 x 10 -11 
(the deep minimum level) and 4 x 10 -11 erg/s/cm 2 (typical 
state), the line followed variations of the power-law compo- 
nent, keeping approximately constant equivalent width. This 
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is clearly a different behaviour than during the ASCA as well 
as the 2001 XMM-Newton observations when, in the same 
range of Fi— 10, the line did not vary, leading to strongly re- 
flection dominated spectra at low-flux states (e.g. Iwasawa 
1996, Fabian & Vaughan 2003, see also figs. 13 and 14 in 
Vaughan & Fabian 2004). This indicates that the nature 
of the physical mechanism underlying variability during the 
2000 observation is different than that of other observations. 
A plausible interpretation involves assumption that, in the 
former, variability was primarily due to variations of intrin- 
sic luminosity while in other cases it was dominated by GR 
effects. 

The rms spectra for the 2001 XMM observation reveal 
the suppressed variability around 6.5 keV, see e.g. figure 16 
in Vaughan & Fabian (2004). Such suppression could be due 
to the presence of reflection from distant material, as sug- 
gested for explanation of similar features in the rms spec- 
tra of several other objects (see, e.g., Ponti et al. 2006). 
In MCG-6-30-15, however, contribution of such a constant 
spectral component is very weak (see Lee et al. 2002, Iwa- 
sawa et al. 1996). Therefore, the suppressed variability of the 
blue peak of the line should be explained within the GR vari- 
ability model. As discussed in Section f3.6l only model Sj<_ T 
(or 5*k) is able to explain such property. Moreover, the de- 
cline around 6.5 keV in the rms spectrum occurs in the range 
of r s for which the model reproduces the time-averaged pro- 
file. 

On the other hand, Vaughan & Fabian (2004) indicate 
hints for increased variability around the Fe line energy in 
the rms spectrum for lower flux states, although ambiguities 
remain in interpretation of that analysis (see their figures 16 
and 17 and discussion in section 5.3.2). We emphasise that 
such property, if confirmed, would contradict predictions of 
model 5k T , where a slightly stronger decline in the rms spec- 
trum should correspond to lower luminosity states, see Fig. 
1151 Then, in the framework of model S^i such property 
would require stronger variations of intrinsic luminosity, at 
a fixed r s , for smaller r s . 

Ponti et al. (2004) note an excess of fractional variabil- 
ity in the 4.7-5.8 keV band (followed by a drop at energies 
close to 6.4 keV) in the rms spectrum of the 2000 obser- 
vation, suggesting presence of a spectral component which 
varies more than the underlying continuum. Such an excess 
may result from variations of the red wing of the line corre- 
sponding to changes of the position within r s — 2 — 4, see 
Fig.ESc). 

4.3 Time-resolved profiles; difference spectra 

High-quality spectra for time-scales on which the flux re- 
mains approximately constant should reveal line profiles 
corresponding to a fixed position of the X-ray source. 
Specifically, for model Sk T these profiles should exhibit 
clear bumpy features in the red wing. MCG-6-30-15 varies 
strongly on time-scales < 1 ks. Current data allow to de- 
termine spectra for time bins of at least 10 ks. Then, these 
spectra would involve averaging over certain ranges of r s in 
model S^ T , which in turn would smooth out the red wing 
(see also Fig. \22{ c)). We note, however, clear indications of 
bumps in the line profile below 6 keV in the ASCA observa- 
tions and some, less significant, hints in the XMM-Newton 
data. If interpreted in terms of model Sk T predictions, this 



would indicate that in some cases generation of X-rays is 
localised around specific r s on timescales > 10 ks. 

The time-averaged profile of 1999 ASCA observation 
(see figure 1 in Shih et al. 2002) has a main peak at ~6.5 
keV and then an additional peak at > 5 keV (qualitatively 
similar but less pronounced peaks around 5 keV are seen in 
time averaged profiles from 1994 and 1997 ASCA observa- 
tions; see Tanaka et al. 1995 and Iwasawa et al. 1999). This 
profile is modelled with a rather flat radial emissivity, q < 3, 
and a small outer radius, in the range r out = 10-30. For such 
parameters, the lower energy peak is formed in the way il- 
lustrated in Fig. 11(c). However, that fitting model requires 
some complex geometry of the central region, with obscura- 
tion of more distant parts of the disc, justifying the above 
values of r out ; otherwise, even a single point-like source in 
the central region would illuminate more distant regions of 
the disc, giving rise to a rather smooth profile, without the 
peak at 5 keV. On the other hand, the observed profiles may 
be explained by increased X-ray activity at r s « 3 in model 
5*k T , without assuming such small values of r out . 

Such clear features are not revealed by XMM-Newton 
for similar flux states. However, the time-resolved (in 10 
ks bins) spectra often show some features at (2-3)<7 level 
around 3.5-4 keV in the 2000 observation (see figure 8 in 
Reynolds et al. 2004; note that the time-average profile in 
that figure has a clear bump with maximum at < 4 keV) 
and around 4.5-5 keV in the 2001 observation (see figure 18 
in Vaughan & Fabian 2004). Finally, Dovciak et al. (2004a) 
note features near E ~ 4.8 keV in the time-averaged line 
profile of 2001 observation, which could be interpreted as a 
Doppler horn resulting from locally enhanced emission on 
some part of the disc. The above energies would correspond 
to r s ~ 2 and 3 for the 2000 and 2001 observation, respec- 
tively. 

Given a rather poor quality of the time-resolved spec- 
tra, less direct approaches to investigate spectral changes 
are made. Specifically, the difference spectra between vari- 
ous flux levels are determined, which show only the variable 
component of the spectrum, while the non- variable compo- 
nent is subtracted away. Figure 15 in Vaughan & Fabian 
(2004) shows the difference spectrum between the average 
spectra, for the whole 2001 observation, with fluxes lower 
and higher than the mean. The difference is consistent with a 
power-law, with residuals from the power-law model not ex- 
ceeding 10 per cent. Therefore, it indicates very weak change 
of the spectrum of the reprocessed component. This prop- 
erty strongly challenges the GR models of variability, which 
in general predict rather significant changes of the line pro- 
file resulting from the change of the source position. Below 
we present example difference spectra for models Sk T an d 
C. 

We note, however, that both the low and high flux spec- 
tra in the analysis in Vaughan & Fabian (2004) involve aver- 
aging over wide ranges of flux states, see their figure 12, and 
the average fluxes for high and low flux spectra differ by a 
factor of 1.5. In GR variability models, such averaging would 
correspond to averaging over large ranges of r s . Taking into 
account that, e.g., model S^ T should involve variations of 
intrinsic luminosity by a factor of < 2, these ranges of r a may 
largely overlap. The most interesting constraint on spectrum 
presumably corresponding to an approximately fixed posi- 
tion is discussed in Section f4. 61 
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Figs. I20f d-f ^ show difference between spectra for r s — 3 
and r s — 1.8 (i.e. between extreme positions considered 
above for the time-averaged profile) in model S^ T and for 
hs = 8 and h s — 4 in model C. The difference between the 
total spectra (Fe line + power-law) is divided by the differ- 
ence of the primary power-laws. Variable red wing gives rise 
to an excess at < 5 keV, not exceeding 10 per cent. The exact 
energy of this excess depends on the value of r a attributed to 
the high state, with higher energies corresponding to larger 
r B . Owing to effects discussed in Section [3.5.2l model S^ T is 
characterised by relatively small variations of the blue peak. 
Deviations of the difference spectrum from a power-law, at 
the peak energy, depend on the angular emissivity law, with 
limb-brightened laws giving rise to smaller deviation. 

In models A and C significantly stronger variations of 
the blue peak occur, therefore, these models predict much 
larger residuals around 6 keV (than model <Sk T ) in the dif- 
ference of spectra between positions corresponding to low 
and high fluxes, see Fig. I20f f). Note also that significant de- 
viations from power-law may occur in the difference between 
two profiles which fit the same observation equally well, see 
the dotted curve in Fig.lSuTf). cf. Miniutti et al. (2003). 

The difference spectrum between the 2001 and 2000 
XMM observations shows an excess at 6.6 keV, see figure 
20 in Vaughan & Fabian 2004, which could be due to a 
change in the flux of the blue peak of the Fe line. Such a 
change of the blue peak, yielding appropriate level of resid- 
uals around 6.6 keV, could result from significant change 
of the source height, as illustrated in Fig. I20f f). with the 
2001 observation corresponding to larger h s . On the other 
hand, this interpretation would not be consistent with other 
properties of the 2001 observation, discussed in this paper, 
which seem to require small h B . Moreover, the line in the dif- 
ference between the 2001 and 2000 observations is narrow, 
while that resulting from the change of h s would be rather 
broad, with significant flux below 6 keV. 

4.4 Equivalent width 

High values of both the equivalent width (EW) of the Fe 
line and the normalisation parameter of Compton reflection 
(R — Q/2n), characterising the X-ray spectrum of MCG- 
6-30-15, are supposed to result from enhancement of disc 
irradiation by GR effects. Slightly different values of EW 
for the 2001 observation result from the two XMM-Newton 
detectors. Specifically, using the same model (with R = 1), 
Fabian et al. (2002) and Fabian & Vaughan (2003) find EW 
= 550 eV and EW = 685 eV using the data from MOS and 
pn camera, respectively. Allowing for higher normalisation 
of reflection, R = 2.2, Fabian et al. (2002) obtain EW = 
410 eV. Using the BeppoSAX data, Ballantyne, Fabian & 
Vaughan (2003) find a high value of reflection parameter, 
R > 2.6, in agreement with the above EW values. 

Note that procedure of computing the EW is particu- 
larly important in GR models. Specifically, for high R, the 
EW with respect to the sum of the power-law and Compton 
reflection, denoted below as EWtot, may be much smaller 
than the EW with respect to the power-law only, denoted 
as EWpl (see also examples of EW values for various pro- 
cedures of computing in Miniutti et al. 2004). 

We have determined both EWtot and EWpl (all val- 
ues given below correspond to # bs = 30°). The latter may 
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Figure 21. The solid curve in panel (a) and (b) shows the 
line EW (with respect to power-law + reflection continuum) 
and the reflection parameter, respectively, as a function of the 
source distance r s , for S^ T with /(/x e m) oc 1 + 2.06/x cm . The 
reflection parameter is determined as R =EWpL/EWp£ ss , see 
text, with EW^ SB = 260 eV. The dotted (blue) curve is for 

/(/Jem) OC log(l + 1/fJ.em) 



be used, by comparing with the classical value (i.e. for the 
model with no relativistic effects), EWp£ ss , to estimate the 
reflection normalisation parameter as R =EWpl/EW p '£ ss . 
The value of EWp£ ss depends on the angular emission law, 
in accordance with changes of the line profile discussed in 
Section [321 E.g., EWp'ff 8 = 170 and 260 eV corresponds 
to I(fi em ) oc log(l + l//i e m) and /((j cm ) oc 1 + 2.06/i em , 
respectively. Neglecting the factor 1.3 (related with elemen- 
tal abundance, see Section 2), we find EW P l ss = 130 eV (in 
agreement with George & Fabian 1991; see their figure 14) 
for /(/i cm ) oc log(l + l//iem). However, dependence of R 
on specific 7(/i e m) is rather weak; differences between limb- 
brightened and limb-darkened cases do not exceed 20 per 
cent. In Fig. 1211 b) we show values of R corresponding to 
limb-darkened 7(/i e m), which is closer to angular dependence 
of Compton reflected component. 

In model S^ T with small r B , the EWtot is extremely 
sensitive to r s , see Fig. [2TT a). For I(fi Bm ) oc 1 + 2.06^i C m, 
the EWtot increases from 480 to 1030 eV with the distance 
decreasing from r s = 3 to 1.8. The superposition of sources 
within r s = 1.8 — 3 (Di, the solid curve in Fig. 1201 a) ) yields 
EWtot = 670 eV. For I(fj, e m) oc log(l + l//i cm ), in the same 
range of r s , EWtot increases from 360 to 800 eV, and the 
superposition yields EWtot = 520 eV. Then, for the limb- 
brightened I(fJ, em ), EWtot is in approximate agreement with 
spectral fits of the XMM-Newton data. The limb-darkened 
I(Hem) yields a slightly higher value, although this conclu- 
sion is affected by uncertain abundance of iron. 

The reflection parameter for the superposition of 
sources within r s = 1.8 — 3 in model S^ T is R — 3.8, which 
value is within the confidence limit from analysis of the Bep- 
poSAX data in Ballantyne et al. (2003). 

The EW is significantly reduced if even a small frac- 
tion of X-rays is emitted at larger distances, as required by 
the time-averaged profile for limb-brightened I(fj, enl )- E.g, 
an additional source at the axis, r s = 10, in model [5£ T ,A] 
(the solid curve in Fig. 1201 b)). emitting 10 per cent of the 
X-ray luminosity, yields EWtot = 380 eV and R<3. 
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4.5 X-ray luminosity 

The primary emission, observed at small ^obs, is strongly 
reduced in model S^ T for the range of distances relevant to 
modelling of the line profile and variability effects. In par- 
ticular, for r s = 3 and 1.8, it is reduced by a factor of 8 and 
70, respectively. Then, the observed hard X-ray flux would 
be reduced by over an order of magnitude due to GR effects. 
The 2-10 keV flux in the typical state of MCG-6-30-15 cor- 
responds to isotropic luminosity Lf-io — 4 x 10 42 erg/s, 
implying the total isotropic luminosity Lx° < 2 x 10 43 erg/s. 
The above reduction factors imply that the intrinsic lumi- 
nosity of the X-ray source is Lx — 2 x 10 44 erg/s ~ 0.5Z/Edd 
(for M = 3 x lO 6 A'/ ) and this in turn indicates M > M E dd, 
making the central engine of MCG-6-30-15 closely similar 
to those of Narrow Line Seyfert 1 galaxies (e.g. Collin & 
Kawaguchi 2004). 

The NIR/optical/UV luminosity (presumably thermal 
emission of an optically thick disc) is estimated to Ldisc > 
2 x 10 43 erg/s by Reynolds et al. (1997). The above lower 
limit corresponds to the minimum reddening allowed for 
MCG-6-30-15. The putative rapid rotation of the black hole 
implies that the thermal disc emission, driven either by di- 
rect dissipation of the accretion power or by irradiation by 
X-rays, should emerge mostly from distances of a few _R g . 
Then, the disc emission would be subject to similar reduc- 
tion as the X-rays and therefore the disc luminosity may be 
comparable to that of the X-ray source (or higher for larger 
values of the reddening). Actually, this scenario, i.e. a more 
luminous disc emission beamed along the equatorial plane, is 
consistent with very high luminosity of MIR/FIR radiation 
coming from reprocessing of central emission, see discussion 
in Reynolds et al. (1997). 




Figure 22. The figure shows the difference between spectra of 
higher flux states, averaged over some ranges of source position, 
and the lowest flux state, for a fixed source position, in model S^_ T 
with I(ficm) oc 1 + 2.06/i cm . The solid (black) curves in panel (a) 
and (b) show superposition of profiles for r s = 1.8 — 3, same as the 
solid curve in Fig. rgoTal. The solid (black) and dotted (blue) curve 
in (c) shows superposition for r B = 2 — 2.2 and 2.5 — 3, respectively. 
The dashed (red) curves show profiles for r s = 1.5 (a) and r B = 1.7 
(b,c). Normalisations of spectra for the superposition and for the 
single source position are adjusted so that they correspond to 
equal times of observation. The solid (black) curves in panels (d-f ) 
and the dotted (blue) curve in (f) show the difference spectra (the 
lower flux spectrum subtracted from the higher flux spectrum) for 
parameters indicated for the top panels. The dotted (blue) curve 
in panel (d) is for the same parameters as the solid one but with 
the spectrum for r s = 1.5 multiplied by a factor of 1.5. 



4.6 Low flux state 

The lowest flux states, with spectra dominated by the reflec- 
tion component, are particularly interesting for studying the 
strongly relativistic regime of space-time. In the GR mod- 
els, these states would correspond to positions of the X-ray 
source in the immediate vicinity of the event horizon. Prop- 
erties of the lowest flux spectra seem to deviate between 
various observations of MCG-6-30-15. 

The line profile of the deep minimum state analysed by 
Iwasawa et al. (1996) is completely dominated by the red 
wing, while the blue peak is strongly suppressed. This could 
be interpreted as the result of generation of primary emission 
at small r s close to the axis (i.e. in the regime of model 
A rather than S), however, some challenge for this model 
results from enhancement of the blue peak by illumination 
of the surrounding disc (i.e. beyond a few R s ) by returning 
radiation. 

Indeed, such a line profile, dominated by the red wing, 
can be modelled with very steep £F e profiles (q > 4) on the 
whole disc, so that contribution from beyond a few R s is 
negligible. As a result, the line is formed in bulk by strongly 
redshifted photons and has a bumpy shape without a blue 
peak. A similar (more extreme) effect is considered as the 
reason of non-detection of a broad line in some spectra pre- 
sumed to be affected by strong GR effects, e.g. NGC 4051 
(see e.g. figure 7 and the related discussion in Miniutti & 
Fabian 2004). In principle, such lines may result from direct 



illumination by a source, with V = 0, at a very small r B 
close to the symmetry axis. However, in any such model, 
additional illumination by the returning radiation would re- 
sult in generation of a rather pronounced blue peak. The 
above effect does not occur if the reflecting (optically thick) 
material is absent within at least 2R g . 

Reynolds et al. (2004) argue that the lack of the blue 
peak is the characteristic property of the low state. However, 
some low flux profiles do have strong blue peak, see below, 
on the other hand, a profile with no blue peak was revealed 
during the flare in the 1997 ASCA observation (Iwasawa 
et al. 1999). Note also that the analysis in Reynolds et al. 
(2004) involves a wide range of flux states including much 
higher than that of the deep minimum. Moreover, their fits 
include a narrow Fe line with EW = 100 eV, which value is 
much higher than these typically found for this object (e.g. 
Iwasawa et al. 1996, Lee et al. 2002). 

Finally, we argue that properties of the low flux state of 
the 2001 observation, analysed by Fabian & Vaughan (2003), 
are consistent with predictions of model . The spectrum 
of that state is strongly dominated by reflection, with very 
large EW = 1.4 keV, similar as for the deep minimum of 
the ASCA observation. On the other hand, the profile has a 
strong blue peak, as implied both by consistency with higher 
flux profiles, see below, and by the fit with q = 3.3 and high 
r out . To account for such high values of EW, model Sk T re- 
quires the X-rays to be emitted from r s < 1.7. For these r s , 
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the observed primary X-ray flux is reduced by a factor of 
> 100, in approximate agreement with normalisation of the 
power-law component for low flux state in figure 5 in Fabian 
& Vaughan (2003). Note that for these r s the strongly irra- 
diated part of the disc, under the source, does not contribute 
to the line profile - due to extreme redshift and reduction 
of this emission. Then, the line can be modelled by a sin- 
gle radial power-law, with q ~ 3, in agreement with the fit 
from Fabian & Vaughan (2003). Such line profile does not 
possess any bumpy features in the red wing and is similar 
to the profile averaged over the range of r s corresponding 
to higher flux states. This property is consistent with those 
derived from analyses of the 2001 observation, as discussed 
below. 

The light curve for the 2001 observation in general ex- 
hibits rapid variability, except for a rather prolonged period 
( > 10 ks) of the lowest flux state, see bin 303:k in figure 2 
from Fabian & Vaughan (2003). Weak changes of the flux 
suggest that this period may correspond to a fixed location 
of the source. Remarkably, spectrum from this period is sim- 
ilar to spectrum averaged over some some range of higher 
fluxes; see figure 3 in Fabian & Vaughan (2003) with up 
to 20 per cent residuals in the difference spectrum. Simi- 
lar property is found in Fabian et al. (2002) for another low 
flux period; Vaughan & Fabian (2004) note that they do not 
find evidence for strong systematic residuals in the difference 
spectra between the lowest flux state and several higher flux 
levels. 

Fig. [22] shows the difference spectra between the flux 
states, averaged over some ranges of r s , and the lowest flux 
state, predicted in model f° r limb-darkened angular 

emissivity - which gives rise to the largest deviations in the 
difference spectra. For r s = 1.7, the residuals are small, ap- 
proximately within the observational constraints. For very 
small r s , model Sk T predicts decrease of the line flux, cf. 
Fig. EIc), therefore, r s = 1.5 the residuals are relatively 
large. However, the shape of the line is similar to the aver- 
aged profile and a slight modification of the model, yielding 
higher line flux, reduces the residuals to a small level, see 
the dotted curve in Fig. !22f dV 

Note that in two Narrow Line Seyfert 1 galaxies, NGC 
4051 and 1H 0707-495, properties of the lowest flux state are 
different than predicted by model S^ T . In these objects, the 
reflected component is correlated with the primary power- 
law at low flux values and becomes more constant at higher 
fluxes (Fabian et al. 2004, Ponti et al. 2006) . Similarly as for 
the 2000 observation of MCG-6-30-15, significant changes of 
intrinsic luminosity would be required to explain those low 
flux states. 

4.7 Discussion 

Model Sk T , with a = 0.998, seems to be the most feasible to 
explain these properties of the X-ray spectrum of MCG-6- 
30-15 which are supposed to result from GR effects. Below 
we summarise predictions of the model which are consistent 
with the observed properties 

(i) Primary emission from sources distributed within r s = 

1.8 — 3 gives rise to the line profile possessing a red wing 
consistent with that revealed in the 2001 observation. 

(ii) Changes of the distance of the source, within r s = 1.8 — 3, 
yield variations of the primary emission flux with constant 



normalisation of reflection component. 

(iii) Change of the distance, within r s = 1.8 — 3, gives rise 
to a decline in the rms spectrum at 6.5 keV. 

(iv) Primary emission from a source located at r s < 1.7 pro- 
duces a strongly reflection dominated spectrum, consistent 
with that observed in the lowest flux states; moreover, the 
line profile for such location is similar to the average profile 
for larger r s . 

(v) Rather significant changes of the line profile occur when 
the distance of a single source changes at r s > 2; the result- 
ing, systematic residuals in the difference spectra are, how- 
ever, relatively small (~ 20 per cent level) and the model 
reproduces reasonably well a power-law like difference spec- 
trum. 

Some details of the model, e.g. modelling of the blue peak 
or the EW value, depend on the angular emission law. 

Note that model Sk, with some modifications, could 
explain a smoothed red wing even without superposition, 
i.e. for a single source at r s = 2 — 3. Moreover, such model 
could account for small changes of the line profile corre- 
sponding to the change of r s . Namely, photons forming the 
variable, bumpy part of the red wing may be depleted due 
to ionization of the disc surface, neglected in this paper. 
Qualitatively, the strongest ionization should occur under 
the source, where the photons forming the redshifted horns 
originate. We emphasise that if future observations confirm 
that the line profile remains strictly unchanged at any flux 
level, applicability of model S would indeed require strong 
reduction of fluorescent emission from the region under the 
source, involving ionization of that region. 

Another possibility to suppress the bumps in the red 
wing involves generation of X-rays at a slightly larger h s , for 
which a more extended area under the source is illuminated 
and the wing does not posses a pronounced feature, compare 
Figs.QHe) and Of). 

Our model S|J T represents most closely a model with 
flares occurring randomly at various distances. In a realis- 
tic scenario, the X-rays would be always generated by flares 
distributed over a range of r s and the observed Fe profile 
would correspond to the superposition of profiles for a single 
r s . Variability effects would then result from changes of the 
distance of the strongest activity. These changes of r s giv- 
ing the dominating contribution should result in changes of 
the line profile described above. We note, however, a caveat 
for the model with flares. Namely, the model involving com- 
pact sources corotating with the disc predicts quasi-periodic 
modulation of the primary emission because of Doppler ef- 
fect as the source circles around the centre (see Zycki & 
Niedzwiecki 2005). This signal would appear on the Keple- 
rian time scale, which is shorter than time scales considered 
in this paper. However, such a signal is not observed in the 
XMM-Newton data from MCG-6-30-15 (Vaughan, Fabian 
& Nandra 2003). This may imply a continuous spacial dis- 
tribution of the hard X-ray source, e.g. an extended corona 
covering the disc surface or a small hot torus replacing the 
disc, in both cases within a few innermost R g (as needed for 
the variability effects). 
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5 DISCUSSION AND CONCLUSIONS 
5.1 Variability models 

We have extended the model, formulated by Miniutti & 
Fabian (2004), relating reduced variability of the reflected 
emission to changing magnitude of relativistic effects as lo- 
cation of the primary X-ray source changes. We find that 
original computations of Miniutti & Fabian (2004) - for a 
model involving a vertically moving source - overestimate 
the reduction effect by assuming a value of the outer radius 
of the disc which is too small for the range of the source 
heights considered in that model. 

On the other hand, we find a significant reduction of the 
variability of reflected emission in a model with a rapidly 
rotating black hole and a source moving radially, low above 
the disc surface. The reduced variability occurs then for the 
innermost range of radial distances, < 4R g . We find also 
that - only in this range of parameters - the GR effects give 
rise to a significant decline around 6.5 keV in rms spectra. 

Note that generation of strong X-ray emission at these 
distances (< 47? g ) is consistent with the condition of a high 
value of a, as - for a rapidly rotating black hole - most ac- 
cretion power is dissipated within a few innermost R g . 



5.2 The black hole spin 

Determination of the value of black hole spin or, even more 
fundamentally, verification of effects predicted by the Kerr 
metric solution of GR equations - remains a major issue of 
black hole astrophysics. In this context, several effects are 
taken into account for X-ray spectroscopy. Strong redshift 
of photons forming the observed red wings, at E < 4 keV, 
is considered as evidence of rapid rotation of a black hole 
(e.g., Brenneman & Reynolds 2006). However, similar red- 
shift can be obtained for a = if emission from r<j < 6 is 
taken into account (Reynolds & Begelman 1997). Then, the 
derived high value of a relies on assumption of no neutral 
iron emission from within the radius of marginal stability. 

Response of the line to increase of primary emission has 
been suggested for future studies. In particular, Reynolds 
et al. (1999) indicate a bump occurring in the line profile 
and progressing to lower energies, with proceeding time, as 
a feature characteristic for a rapidly rotating black hole. 
Again, a similar - redshifted and Shapiro delayed - bump 
should appear for a non-rotating black hole if fluorescence 
inside r ms was taken into account. 

A straightforward analysis of the space-time metric 
could be performed for systems observed close to edge-on, 
for which effects due to lensing by a black hole would be di- 
rectly seen in the line profile (e.g., Zakharov & Repin 2003). 
However, as noted in Narayan & McClintock (2005), there 
seems to be a selection effect preventing such systems from 
being observed. 

Then, we point out that a (largely unambiguous) anal- 
ysis of imprints of the space-time metric would be possible 
in profile of the line resulting from irradiation by a strong 
flare just above the disc surface. If such a flare occurred at 
r s < 4, properties of the space-time related with the black 
hole rotation would result in > 10 per cent in magnitude 
effects in the line profile. The effects related to the value 
of a are rather subtle but in principle possible to establish 



observationally. A compact flare dominating total emission 
would be required to make such an analysis feasible and vi- 
ability of such scenario is uncertain. Such flares are indeed 
occasionally observed (e.g., Ponti et al. 2004). Interestingly, 
however, the Fe Ka line was actually very weak during the 
flare analysed by Ponti et al. (2004), and a strong line ap- 
peared in the spectrum with significant time delay after the 
flare. 

Finally, we emphasise that the reduced variability of 
reflected component - basing on mechanism advocated in 
this paper - is itself a direct manifestation of the nature of 
the Kerr space-time. Note that another effect resulting from 
properties of the Kerr metric, namely anisotropic emission 
of hard X-rays generated close to a rotating black hole, is 
qualitatively consistent with inclination-angle dependence 
of intrinsic spectra of Seyfert galaxies (Niedzwiecki 2005). 
In support for the tentative relation of these two effects to 
properties of the Kerr metric, note also that comparison of 
the total mass in black holes in the local universe with the 
total luminosity produced by active galactic nuclei indicates 
that most supermassive black holes should rotate rapidly, 
e.g. Elvis, Risaliti & Zamorani (2002). 

5.3 Modelling the red wing of relativistic Fe lines 

As discussed above, detection of photons with E < 4 keV 
is considered as the evidence of rapid rotation of the black 
hole. Such strongly redshifted photons were revealed in the 
Fe line profiles observed in MCG-6-30-15 and several other 
objects (e.g., Miller et al. 2004, Miniutti et al. 2004). When 
fitted by models assuming a power-law radial emissivity, the 
observed profiles require q > 4 within (6 — 10) R g . Consider- 
ing physical scenarios for generation of such profiles, we find 
that they 

(i) can be produced as a result of illumination by a hard 
X-ray source located close to a black hole (preferably at 
r B — 2-3) and rather close to the disc surface (h B < 1); 

(ii) cannot be explained by models involving a source at a 
height of several R g , or higher, close to the axis; 

if they are to come from a disc, with small inclination, sur- 
rounding a rotating black hole. Regarding (ii) we find, in 
agreement with previous studies, that such location of the 
source yields a steep £f c , with q > 3, in the innermost part 
of the disc. However, this steep ep c occurs only within < 2 
(and due to gravitational blueshift rather than light bend- 
ing). The majority of emission from this region of the disc is 
either captured or bent toward the disc plane; moreover, the 
observed photons are strongly redshifted. As a result, this 
emission gives only a minor contribution to the line profile 
and it is not relevant to the shape and strength of the ob- 
served red wings. This conclusion is not affected by effects 
related with azimuthal motion of the source, returning ra- 
diation or angular emission law of fluorescence; these effects 
influence mostly the strength of the blue peak. On the other 
hand, the red wing is primarily related with position of the 
X-ray source. 

5.4 Returning radiation, non-local irradiation 

A further challenge for modelling certain line profiles re- 
sults from illumination of the surrounding disc (i.e. beyond 
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a few 7?g) by returning radiation. We find that this effect 
may be strong for the range of parameters (small r s and V) 
previously not explored. A similarly strong enhancement of 
Compton reflection by returning radiation, for very small 
distances (r 8 < 2) of primary source, was noted recently by 
Suebsewong et al. (2006). 

Another effect, related to bending of photon trajecto- 
ries to the disc plane, affects models relating the Fe radial 
profile to some physical processes. These models typically 
assume that the radial profile of Fe emission is the same 
as some physically motivated profile of energy dissipation. 
E.g., Reynolds et al. (2004) make such assumption in their 
analysis of the XMM observation of MCG-6-30-15, applying 
model of the torqued-disc emission (where emission results 
mostly from extraction of rotational energy of the black hole 
and thus is very centrally concentrated). Moreover, they as- 
sume that there is no dissipation, and thus no reprocessing, 
beyond rather small (several R g ) outer radius. We empha- 
sise that even for X-rays generated very low above the disc 
surface, a significant illumination of more distant regions of 
the disc should occur. 

Both the returning radiation and the non-locality of ir- 
radiation result primarily in enhancement of the blue peak 
of the line. 

5.5 Azimuthal motion 

Impact of strong gravity is usually studied under specific 
assumptions on the X-ray source motion. Typically, angu- 
lar velocity of the source is related to that of the disc (e.g. 
Ruszkowski 2000; Miniutti & Fabian 2004), while Dabrowski 
& Lasenby (2001) assume a static primary source. We note 
that change of V significantly affects flux and profile of 
the iron line through combination of SR and GR effects. 
Moreover, the Kerr metric terms yield a non-trivial depen- 
dence between V and SI (equation (©)■ Then, the assumed 
parametrisation of motion may be crucial for derived prop- 
erties, e.g. in some variability models. Obviously, additional 
effects may result from relativistic vertical or radial motion 
of the source (e.g. Yu & Lu 2001), which were neglected in 
this paper. 

5.6 Applicability of our results 

We focused above on MCG-6-30-15, for which several ob- 
served properties may be explained by our model. Simi- 
lar effects, including reduced variability of reflection, pro- 
nounced red wings or strongly reflection dominated spectra, 
have been revealed in some Narrow Line Seyfert 1 galaxies 
(e.g. Fabian et al. 2004, Ponti et al. 2006), as well as in stel- 
lar mass black-hole systems (e.g. Miniutti et al. 2004) and in 
galactic nuclei with much higher black hole masses (Fabian 
et al. 2005). Note that our scenario implies that the X-ray 
luminosity is underestimated by an order of magnitude, in 
these objects. 

On the other hand, most black-hole systems do not show 
signatures of such extreme effects. Then, similar reduction 
of the X-ray flux does not necessarily characterise any low 
inclination system. Note, however that such property would 
be inevitable for rapidly rotating black holes, where major 
fraction of accretion power is dissipated at very small r s , if 
this power is converted into hard X-rays in situ. 



As in virtually all previous studies of that subject, we 
analysed impact of GR effects for the simplest case of an 
isotropic point source of hard X-ray emission, which ap- 
proach follows from our lack of understanding of the X-ray 
source nature. The derived properties result from transfer 
of radiation from source to disc and observer and from disc 
to observer. Additional (but smaller in magnitude) effects 
may occur in more realistic models. E.g., Comptonization 
close to a Kerr black hole gives rise to anisotropic emis- 
sion, cf. Niedzwiecki (2005), therefore radiation with differ- 
ent spectral index may irradiate various parts of the disc, 
while the transfer effects affect only normalisation and cut- 
off energy of the primary emission. 

The simplified description of X-ray source considered in 
this paper approximates most closely scenario with magnetic 
flares above the disc surface. Then, our results may be di- 
rectly applicable to a model with flares occurring randomly 
at various radial distances. Qualitatively, we may assess sim- 
ilar variability effects for continuous spacial distributions of 
the hard X-ray source, e.g. an extended corona covering the 
disc surface or a small hot torus replacing the disc within 
a few innermost R s . Varying size of such an extended, hot 
plasma, in the Kerr metric, should give rise reduced variabil- 
ity of reflected component. In particular, a decline of direct 
emission from a shrinking corona or torus would be observed, 
even if its intrinsic luminosity remained unchanged, while in- 
creasing fraction of its emission would be bent to the disc 
plane giving rise to strong reflection component. 

The major shortcoming in our study results from the 
neglect of ionization effects. The strongest ionization should 
occur below the source, especially in models with low height 
above the disc surface. In general, strong ionization of this 
region should suppress the redshifted Doppler horns formed 
in the red-wing. On the one hand, this would make studies of 
effects related to value of a less feasible. On the other hand, 
such depletion of the variable contribution to the red wing 
would reduce variations of the line profile at various flux 
states. However, details of ionization structure, and of the 
related reduction of the variable contribution to the lien, de- 
pend on additional assumptions, in particular, on azimuthal 
distribution of flares. Obviously, a strong single flare would 
give rise to stronger ionization than many weak flares uni- 
formly distributed at a given r s . 



5.7 Summary 

If attributed to strong-gravity effects, both the time- 
averaged line profile and the variability pattern observed 
in MCG-6-30-15 by XMM-Newton independently indicate 
that a primary hard X-ray source must be located very close 
(<4_R g ) to a black hole, i.e. in the region where Kerr met- 
ric effects become crucial. Rapid rotation of the black hole 
is necessary (and vastly sufficient) to account for reduction 
of variability of reflected emission for such spacial location 
of the source. Bending to the equatorial plane, underlying 
this reduction effect, appears to be the most pronounced ef- 
fect of the Kerr metric to be studied in the X-ray spectra of 
black-hole systems. 
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